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On the shapes of elementary domains or 
why Mandelbrot Set is made from almost ideal circles? 

V.Dolotin and A.Morozov 
ITEP, Moscow, Russia 
ABSTRACT 

Direct look at the celebrated "chaotic" Mandelbrot Set in Fig. Q] immediately reveals that it is a collection of 
almost ideal circles and cardioids, unified in a specific forest structure. In /hep-th/9501235 a systematic algebro- 
geometric approach was developed to the study of generic Mandelbrot sets, but emergency of nearly ideal circles 
in the special case of the family x 2 + c was not fully explained. In the present paper the shape of the elementary 
constituents of Mandelbrot Set is explicitly calculated, and difference between the shapes of root and descendant 
domains (cardioids and circles respectively) is explained. Such qualitative difference persists for all other Man- 
delbrot sets: descendant domains always have one less cusp than the root ones. Details of the phase transition 
\ between different Mandelbrot sets are explicitly demonstrated, including overlaps between elementary domains and 
£NJ 1 dynamics of attraction/repulsion regions. Explicit examples of 3-dimensional sections of Universal Mandelbrot Set 
£H , are given. Also a systematic small-size approximation is developed for evaluation of various Feigenbaum indices. 
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1 Introduction 

The question of how dynamics of a physical system depends on the choice of its Hamiltonian is one of the most 
important in theoretical and mathematical physics. Its significance is only enhanced by the fact that in modern 
theory dynamics is considered not only in physical time, but in many other variables, including the coupling 
constants of a theory and the shape of functional-integration domain (the so-called renormalization-group dynamics 
[I]). Normally dynamics is described in terms of a phase portrait or of eigenstates configuration for classical and 
quantum systems respectively, and the question is how these portraits and configurations change under variation 
of the Hamiltonian. It is well known that this change is not everywhere smooth: at particular "critical" or 
" bifurcation" points in the space of Hamiltonians the phase portraits get reshuffled and change qualitatively, not 
only quantitatively: this phenomenon is also known as "phase transition" . Normally these bifurcations are described 
in terms of the change of stability properties of various periodic orbits (including fixed points, cycles and "strange 
attractors"). More delicate information is provided by the study of intersections of unstable orbits, but it is a little 
more difficult to extract. 

Dynamical systems are much better studied in the case of discrete dynamics: this reveals many properties, 
which get hidden in transition to continuous evolution. In other words, this resolves ambiguities of continuous 
dynamics: there are many different discrete dynamics behind a single continuous one, and to reveal the properties 
of the latter it is often needed to look at the whole variety of the former. In classical case discrete dynamics (with 
a time-independent "Hamiltonian") is the theory of iterated maps: 

x - m - r\x) = f(f(x)) -...-> r( X ) = /(r^Hx)) - . . . (i.i) 

and is actually a branch of algebraic geometry [5] (for generalization of [5] to discrete dynamics of many variables 
see sections 7 and 8 of [3]). According to [2J, the structure of the phase portrait is controlled by the Julia set: 
collection of all periodic orbits of the map / in the x space, i.e. of all roots of all functions 

F p (x)=f°P{x)-x (1.2) 
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Therefore the Universal Mandelbrot set (UMS), consisting of all points f(x) in the space of functions (Hamiltonians) 
where some two periodic orbits coincide, can be alternatively characterized as the Universal Discriminantal variety 
formed by the roots of various resultants R x (yF pm (x) / F p {x), F p {x)\ . This is almost a tautological identification, 

since by definition the resultant of two functions vanishes whenever they have a common zero, still it establishes 
relation between a priori different sciences: the theory of phase transitions and algebraic geometry. 

Usually considerations are restricted to particular sections of the Universal Mandelbrot set, by choosing specific 
one-dimensional families of functions: see Figs[T][3]for the three famous examples, f(x; c) = x d + c with d = 2, d = 3 
and d = 4. In Fig|4]we show also the result of a deviation from this simple form. We keep the name "Mandelbrot 
set" for any of such one-complex-dimensional sections of infinite-dimensional UMS, while "Mandelbrot Set" (with 
two capital letters) refers to the original example in FigJTJ Today all kinds of experimental data about these sets 
can be obtained with the help of available computer programs, like Fractal Explorer which is used to make 
Figures [T][3] in the present papeiQ. 



Mandelbrot sets are often considered as typical examples of " fractal structures" , serving mostly for admiration, 
philosophical speculations and, perhaps, numerical exercises. However, as explained in [2J, they can actually be 
subjected to systematic scientific investigation, in the style of experimental mathematics, with questions coming 
from direct observations and numerical experiments, and rigorous answers provided by knowledge of underlying 
algebro-geometric structures. Our presentation below can be considered as an example of this increasingly impor- 
tant approach to modern mathematical physics problems. 

As explained in [2] - and clearly seen in FigsQJH - the Mandelbrot set consists of infinitely many separated 
clusters, of which only the central one is well seen in the main picture, while examples of smaller clusters are shown 
in auxiliary pictures with enhanced resolution. Though separated, clusters form a well organized structure: they 
are connected by "trails", populated with other clusters. Further, each cluster has its own tree structure, FigfT]D, 
with two types of elementary domains: one type at the root of the tree and another type at all higher nodes (we 
call them descendants). FigH] demonstrates that a given Mandelbrot set can contain different types of clusters, 
while for special families f(x;c) = x d + c, where maps possess additional Z^_i symmetry x — > e 27Tln /( d ^ 1 ) Xj a n 
clusters are of the same type. In Mandelbrot Set (i.e. for d = 2) the root elementary domains are nearly ideal 
cardioids, FigOA, 

c-c =r e i4, ( 1- ] (1.3) 



while descendants are nearly ideal circles 

c-co = r e i4, (1.4) 
For d > 2 we have nearly ideal (d — l)-cusp cardioids, Fig [51 

c-c = r a e l * ^1- 2 ei(d ' 1) J ( L5 ) 

at roots, while descendant domains are deformed cardioids with d — 2 cusps and some non-vanishing coefficients 
in 

c-c = roe* (l + £ a k e ik ^ - ^ S d ~^*\ (1.6) 

Actual shapes of the domains slightly deviate from these ideal cardioids and circles and depend on particular cluster 
and node, but deviations are at the level of a few percents at most. 
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Figure 1: A. Mandelbrot set M2 for the family f(x\c) = x 2 + c [5]. The boundary of the black domain in the complex plane 
of c-variable consists of all values of c where Julia set is reshuffled: as explained in [2] this happens when a stable orbit crosses an 
unstable one. A given orbit O is stable within and elementary domain, which is - with incredibly good accuracy - either a cardioid 
c— c<p = roe'* (l — ^e 1 *) at the center (root) of a cluster or a circle c— cq = TQe 1 ^ at the non-root nodes of the tree. Explanation of 
this fact (that only these two shapes occur and exactly in these roles: at roots and higher nodes respectively) is the task of the present 
paper. Projection to the line of real c is also shown. B. Two pieces of the Mandelbrot set M2 under microscope: exactly the same 
structures are seen as the central cluster in A, with the same central cardioids and attached circles. There are infinitely many such 
structures of different sizes r<p in M2, since r<p are very small, they are not actually seen in A, but can be easily studied with the 
help of the Fractal Explorer [4]. More numerical characteristics of the lowest elementary domains are collected in a Table in s|6] C. 
Domains (1) and (2, 1), obtained as numerical solutions of exact equations (12.11 1. see s !3.3.2l They coincide with domains, seen in A. 
This is non-trivial, because pictures A and B are obtained by absolutely different procedure (actually, black region consists of points 
c, with limited sequences /® n (c)), and there is no a priori reason for any parts of them to satisfy any kind of algebraic equations. 
Separation of Mandelbrot sets into domains, possessing an algebraic description, in particular the relation between pictures A and 
C, is important property of iterated maps. D. Tree structure of the central cluster (only a few lowest branches are shown), each 
branching occurs at the center of a new elementary domain, and the number at the vertex is the order of the periodic orbit, which is 
stable inside this domain. Thus elementary domains are naturally labeled by sequences of divisors, leading to the root of the tree. All 
other clusters are represented by exactly the same trees, only numbers are multiplied by the order of the root orbit. Thus entire M2 
has a natural forest structure. 
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Figure 2: Mandelbrot set Mz for the family f(x; c) = x 3 + c. Everything said about M2 is true in this case, only the place of simple 
cardioids at roots of clusters is taken by the two-cusp ones c — c<p ss rpc 1 * (l — |e 2l *j . Descendant domains are nicely approximated 
by single-cusp cardioids c — cq rpe 1 ^ (l — ^e l( ^ , see eg . J10.5I I . No circles are present. 




P] l» 

Figure 3: Mandelbrot set A^4 for the family /(x; c) = x 4 + c. Everything said about A-f 2 is true in this case, only the place of simple 
cardioids at roots of clusters is taken by the three-cusp ones c — cq k, rge'* — ie 3l< ^). Descendant domains have the shape of 
deformed 2-cusp cardioid, c - c = r Q e'* (l - ae^ + |e 2,; ^) with a = 2" 4 / 3 « 0.40. . . and 6 « 11/(9 ■ 2 2 / 3 ) = 0.77. . ., see eq. fTuTTt . 
No circles or single-cusp cardioids are present. 
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Figure 4: Mandelbrot set M3—1 for the family f(x;c) = ax 3 + (1 — a)x 2 + c with two different values of additional parameter 
a = 4/5 (A) and a = 2/3 (B). Everything said about M2 is true in this case, only in addition to the simple single-cusp cardioids 
c — co = roe 1 '*' — ^e"^) in the role of central (root) domains there are also two-cusp ones, c — cq = roe 1 ^ (l — ^ e2!< ^) ■ Moreover, 
for distinguished value o = 1, see Figf2] simple cardioids do not appear as central (root) domains of individual clusters at all - their 
place at roots is taken by 2-cusp curves. Instead for a = 1 the single-cusp cardioids fully replace circles in the role of descendant 
domains. For smaller values of a the simple cardioids start to emerge as roots (and circles - as descendants), but in remote clusters, 
at large distances from the central domain. The central cluster in these pictures look a little asymmetric: this is wrong, and is an 
artefact of the erroneous algorithm, used to construct Mandelbrot sets by Fractal Explorer, see introductory remarks to sf4] below. 
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Figure 5: Cardioid curves described by the equations l|l,3|l and ((3J. These pictures reproduce the shapes of the central domains in 
Mandelbrot sets M d . A. d = 2. This is the central domain of M2 in Figfl] B. d = 3. This is the central domain of Mz in Fig[2] 
C. d = 4. This is the central domain of M4 in Fig|3] 



Each elementary domain is associated with some periodic orbit O of the map /, which is stable exactly within 
the domain. Thus the domain is labeled by the order p of this orbit. This fact allows one to give an explicit 
analytical description of the domain's shape, see eq. (|2.1|) below. However, there are many different periodic orbits 
with the same p and thus many elementary domains with the same p. They differ by the choice of the root c p 
of the equation f p (c) = f op (x = x cr ;c) = x cr , which lies in the "center" of the domain (x cr is the root of f'{x), 
f'(x cr , c) — 0, in most examples below x cr = 0). Some of these order-p domains are at roots of some clusters, some 
at higher-level nodes, but in that case the root of the cluster should be associated with a divisor of p. Actually the 
tree underlying the cluster is the divisor tree, and the entire forest structure (i.e. collection of trees, associated with 
all clusters) is that of divisor forest of natural numbers. Accordingly, elementary domain is labeled by a sequence 

of integers (^(p r ■ mi ■ . . . ■ mk) ar k , . . . , (p r ■ m\) ar 1 , (Pr)a r j , to be read from rights to left. Sequence of multiples 

[mfc, . . . , mi] characterizes domain's position in the tree, k is the distance from the root, p — (p r ■ mi ■ . . . ■ mk) is 
the period of the corresponding orbit, andp r is the order of the orbit, associated with the root domain. Since there 
can be many different root domains with the same p r in particular Mandelbrot set and many different descendants 
with the same p in a given tree, there are additional labels a, distinguishing between different trees in the forest 
with the same order p r at roots and between different branches with the same p in each tree. While divisor trees 
are the same for all particular Mandelbrot sets, collections of a's arc different: they are defined by the way the 
given section crosses the p- variety in the Universal Mandelbrot set, since it can be crossed many times, there are 
many traces of the same variety in the given section (in the role of either root or descendant domains) - and this 
is the origin of the forest and of the a parameters, which at the level of particular Mandelbrot set look somewhat 
arbitrary. At the level of UMS there is a single divisor tree and a section intersects it many times and cuts it into 
many similar trees: any cut-off branch looks like a separate tree and gives rise to a separate cluster. The memory 
of their common origin at UMS level is preserved in the trail structure, connecting the clusters, but its detailed 
description is not yet available. 

Two adjacent elementary domains touch at exactly one point (i.e. along a complex-codimension-one variety in 
UMS), where their corresponding orbits cross and exchange stability. This important statement, however, needs to 
be treated carefully: as we shall see, in general (beyond the x d + c families) the elementary domains can overlap: 
there can be several stable orbits at the same value of c. This means that arbitrarily chosen c is not a good 
coordinate on a Mandelbrot set, which is actually a fibration over a region in the complex-c plane than a region 
per se. Fibration structure is inherited from Julia sheaf over the Mandelbrot set [2]. In this general situation 
the word touch is not fully informative: when different domains seem to overlap, they rather lie in different fibers 
over the same region on c-plane, and these fibers are sewed exactly at a single point, where the orbis cross. What 
we can show in simple 2d pictures are projections of the domains, these projections can overlap and touch at the 
orbit-crossing point. Touch means that the tangents to two domains are collinear, in practice this can be either a 
smooth touching (typical for crossing of orbits of different orders) or a cusp (when the orders are the same). 

Crossing of orbits is possible only when the order of the smaller one (with domain lying closer to the root in the 
cluster) divides the order of the bigger one. Analytically, crossing takes place at the root of associated resultant. 
If the two orders differ by a factor of 2, this is the celebrated period-doubling bifurcation [H [8] - and the chain 
of exactly such bifurcations occurs along the real line in Fig[TJ- but in fact doubling is in no way distinguished: 
bifurcation can multiply period by arbitrary integer m. Crossing of unstable orbits is not seen at the level of 

1 However, one should be careful in using this program for non-canonical families, like in Fig|4] see [2] and s|4]below for explanations. 
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Mandelbrot sets shown in FigsQ|3l- to study these phenomena (also essential for bifurcations of Julia sets) the full 
(or Grand) Mandelbrot set should be considered. Actually, behind UMS there are even more fundamental entities: 
the Universal Julia Sheaf (UJS), consisting of all periodic orbits of all orders "hanged" over the UMS, and the 
Grand UJS, including also all pre-orbits of periodic orbits. UMS is a projection of UJS, obtained by neglect of 
the phase-space dimension x, where the orbits live, and, as any projection, it can and does suffer from overlap 
ambiguities, namely, when two different stable orbits coexist at the same point of UMS. 

We refer to [3j for further details and explanations. The task of this paper is to provide close-to-Earth illus- 
trations of somewhat abstract formulations from the previous paragraphs and to fill some of the gaps left in [5] , 
which concern three closely related subjects: 

(i) the shape of elementary domains; 

(ii) Feigenbaum indices, defining the ratio of sizes of adjacent elementary domains (immediate descendant as 
compared to its parent) from the ratio of the corresponding orders; 

(iii) reshuffling of Mandelbrot set and its elementary sets under the change of selected family of functions, i.e. 
new properties of 2c-dimensional sections of Universal Mandelbrot set as compared to the lc-dimensional sections. 

In fact these subjects capture the main aspects of the general theory and at the same time they can be considered 
by simple methods of theoretical physics with minimal involvement of abstract algebro-geometric constructions. 
Even resultant and discriminant analysis, which was the main machinery in [2], will be at periphery of our simplified 
presentation in this paper. 

As a key puzzle and a starting point for all considerations we choose the question, posed in the title of this 
paper: why exactly the cardioids (|1.3[) and circles (|1 .4[) and exactly in the right places in the divisor forest emerge 
as the shapes of elementary domains of the Mandelbrot Set in Fig[TJ and how this picture is continuously deformed 
into Figsl5]and[3 

As already shown in [2], cardioids (|1.5|) exactly describe the central domains (1) for the families x d + c, while in 
description of all other domains they straightforwardly appear in the "small-size approximation" (SSA) to exact 
shape-defining eqs. (|2.ip . In what follows we 

- explain, why for the x d + c families the cardioids (|1.5p do not exhaust all possible shapes: deformed cardioids 
(|1.6j) with one less cusp are also allowed; 

- explain, why ljl.5|l appear exactly at roots of clusters, while all descendants have one cusp less: instead of this 
lacking cusp a descendant domain has a merging-point to the parent domain; 

- provide a detailed description of interpolation between Mandelbrot sets and Julia sheafs for the families x 2 + c 
and x 3 + c (actually only the orbits of two lowest orders p = 1 and p = 2 will be analyzed, but this is already 
enough to reveal many interesting details of the process); 

- demonstrate inaccuracy of Fractal Explorer [4] (and thus the underlying text-book interpretations of Man- 
delbrot sets) in application to UMS studies and appeal for the writing of corrected and fully adequate computer 
programs on the base of improved knowledge provided by [2]; 

- demonstrate high accuracy of the small-size approximation (SSA) in evaluation of various characteristics of 
the Mandelbrot Set by comparing its predictions for the (complex-valued) sizes tq and Feigenbaum indices with 
exact answers (when they are already available) and experimental data provided by the Fractal Explorer [4] . 

Concerning the last story - about the SSA - it deserves saying that no theoretical explanation for its spectacular 
accuracy is known: particular corrections are not small, but various corrections always combine into a small 
quantity, whenever characteristics of Mandelbrot Set are evaluated. At the same time, as explained in [2], SSA fails 
completely in description of Julia sets; still it describes reasonably well the Mandelbrot sets for the families x d + c 
with d > 2 (though some qualitative properties are spoiled, e.g. cusps are somewhat smoothed), but works much 
worse for interpolations between different d. In any case, today SSA remains the only available tool for theoretical 
investigation of high branches in divisor tree, in particular for approximate evaluation (rather than measuring) 
of various Feigenbaum indices - and for this purpose it works spectacularly well, even for d > 2 and even for 
interpolations. Still rigorous algebro-geometric theory of Feigenbaum indices remains to be found. 

We begin in s [2] from qualitative description of elementary domains supported by the limited amount of exactly- 
solvable examples in s|3j these include some non-trivial cases and, as usual, provide the solid ground for future 
approximate considerations. Section [4] is devoted to interpolation {ax 3 + bx 2 + c} between the two best-known 
Mandelbrot sets: {x 2 + c} and {x 3 + c}. Other examples of 3ij-dimensional sections of UMS (actually of its central 
domain) are also given in this section. Then in s[5]we introduce the small-size approximation and present some 
calculations for the Mandelbrot Set. Their results are compared with experimental data in sj6j After some more 
borrowing SSA calculations in s|7j we discuss Feigenbaum indices in s|8] Section[lO]is devoted to SSA consideration 
of some other Mandelbrot sets. Brief conclusions are collected in s fTTI 



8 



2 The shape of elementary domains. Qualitative description 



2.1 Defining equations 

According to eqs.(10) and (38) of [2], the boundary of an elementary domain satisfies a pair of equations: 

Here F p (x) = f op (x) — x, prime denotes x-derivative and 9 is a real-valued angle-parameter used to coordinatize 
the boundary of the domain (it can actually vary between and a multiple of 2n, see below). After x is excluded 
from the pair of equations (|2.ip . we obtain a rcal-codimension-one hypersurface in the space of functions /, i.e. a 
collection of curves c(9) in l^-dimensional Mandelbrot set. Particular curves - branches of c{9) - are boundaries 
of particular elementary domains of the order p, root and descendant. 



2.2 Cusps 

Even if function c{9) is smooth, the corresponding curve in the complex-c plane can be singular. Generical 
singularity is self-intersection, which takes place when c{6\) = 0(62) for 9\ ^ 62- Of interest for us are cusps: 
degenerated self-intersections, appearing in the limit when 82 — ► 9%, i.e. when ^§(#o) = at some #o- In the 
vicinity of such point <r($) = c(6) — c(9q) = ad 2 + bd 3 + . . ., where i9 = 9 — 9q. This means that 



Re 



(£) = tf 2 + . . . , while Im (^) = Im \± J $ 3 + . . . (2.2) 
i.e. 

-©—©war P.3) 

Thus we see that a cusp emerges at points where dc/d9 = and its orientation in the complex-c plane is defined 
by the phase of the complex-valued parameter a. 
If Im(6/a) = 0, i.e. 

(^) 



dc/dO 

along with dc/d9 — 0, then a self-intersection point collides with the cusp and disappears. 



2.3 Cardioids 

Cardioids are represented by polynomials of the unimodular variable, they form the simplest natural class of curves 
with cusps. 

For quadratic cardioid, 

c = r(e l * + ae 2 ^) = £■ ((1 + 2ae^) 2 - l) , (2.5) 

derivative vanishes, dc/d(f> — 0, when 2ae % * = —1. This never happens if \a\ 7^ \. Thus a cusp (and exactly one) 
occurs only when \a\ = i, the curve is everywhere smooth for \a\ < | and possesses one self- intersection for \a\ > i. 



For cubic cardioid, 

c = r (e** + ae 2 ^ + be 31 *) , (2.6) 

derivative dc/ d<$> = vanishes when 

1 + 2ae i<t> + 36e 2i * = 0, (2.7) 

i.e. when the r.h.s. e~ % * — —a ± ^/ a 2 — 36 has unit modulus. If a and b are real, then cusp can occur when either 
1 + 2a + 36 = (then there is one cusp at (f> = 0, unless a = and 6 = — ^, when another cusp appears at (f> — it, see 
FigO or — 1 < a < 1, 6=+l/3 (then two cusps arise at <f> = ±0o 7^ 0, 7r). In general, (|2.7p defines a hypersurface 
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of real codimension one in the space of complex parameters a and b (parameterized by </>), where cubic cardioids 
(|2.6jl have cusps (one or two). 

Transition point (|2.4p between a phase with and without self-intersection is defined by a system of two equations, 

( dc/dcf) ~ 1 + 2a + 3b = 
\ d 3 c/d</> 3 - 1 + 8a + 27b = 

i.e. a = —4/5 and 6=1/5, see Fig|6l 

MAPLE program for cardioid studies, which was used to generate Figsj5] and El can be found in Appendix to 
this paper, see s ll2.ll 



2.4 Cusps in the boundaries of elementary domains 

The second component of eq. ()2.ip implies that 



pdc__ F ,dx 

p de~ p d0 



(dot and prime denote c- and x-derivatives respectively), so that 35 = when F p = 0, provided F p ^ at the 
same point. Together with the first eq. (|2.1[) this means that cusp can occur only when 9 — 0. Thus the number of 
cusps depends essentially on the range of variation of 0- variable. If 9 runs from to 2i:(d — 1), we can expect up 
to d — 1 cusps to occur. 



2.5 Why descendant domains have one cusp less than the root ones? 

Descendant domain differs from the root one, because it always has one special point at the boundary where F' = 
and F = together. This means that there is no cusp at this point, and if the total number of zeroes of F' at the 
boundary was d — 1, but the domain was a descendant, then the total number of cusps will be d — 2. 

Characteristic feature of any descendant is reducibility of the corresponding function F mp : it is divisible by F p 
of a parent domain, 

Fmpi'E*) — GmpfxjFp^), (2.9) 

(in variance with G mp from ref.[2] such G mp can still be reducible, but this does not matter for our purposes in this 
paper). Then F' — G' mp F p + G mp F p and F mp = G mp F p + G mp F p vanish simultaneously whenever both F p = 
and G rnp = 0, i.e. when x belongs simultaneously to orbits of orders p and mp. According to [5] the last two 
equations possess exactly one common zero at the boundary of descendant domain: it is exactly the merging point, 
where descendant domain is attached to the parent one, and in the c-space it is a zero of the resultant R(G m p, F p ). 

Discriminant D(G mp ) also vanishes when R(G mp , F p ) = 0, because different points of the order-mp orbit (roots 
of G mp ) should merge m-wise to merge with the points of the order p-orbit (roots of F p ). Of more interest are 
other zeroes of D(G mp ), representing crossings of different orbits of order mp. 



3 Exact results about elementary domains 

This section is devoted to exactly-solvable examples. Here exact solvability means not obligatory explicit analytical 
solutions - though they will also be considered. Whenever the problem can be effectively studied by user-friendly 
computer tools like MAPLE or Mathematica, it is considered equally (and may be even better) solvable as if explicit 
formulas were derived. We shall see that sometime the best way to analyze such explicit formula is to generate its 
plot with the help of the same MAPLE. One should keep in mind, however, that the number of problems solvable 
in this way is also very limited: in most cases even clearly formulated algebraic problems can be handled only by 
specially designed programs, which usually could but never were written. This makes such problems potentially 
solvable (as many other hot problems in theoretical physics), but they are clearly different from practically solvable. 
We also distinguish these solvable problems from those which are effectively solved, but only approximately: under 
certain additional assumptions or when improving of accuracy is increasingly difficult (like it happens, for example, 
in perturbation theory). We turn to approximate methods in sslSTITOl Before we are going to describe what is 
known today at exact level. 
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a=L/5 t b^X'IS 




a=-l/8. b=-l/4 a =G, b=-j/3 

Figure 6: Cardioids with cusps and self-intersections. The cusp-possessing subset in the family 1 12,61 1 with 1 + 2a + 36 = is shown. 
A. a = i, b = — j£. Both cusp and self-intersection are present. B. a = — A, b = i. This is the point where self-intersection point 
hits the cusp and disappears. C. a = — ~, b = — j. A single cusp is present at </> = 0. D. a = 0, b = — ~. The second cusp appears 
at (p = TV. 
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Our primary goal is to understand what is the domain of variation of the 9- variable - because we already know 
from s !2.4l that it is its size that defines the number of cusps, both for root and descendant domains. Moreover, 
we want to see how this variation domain is changed in transition from one Mandelbrot set to another, i.e. to 
study the bifurcations of Mandelbrot sets themselves, which happen in complex-codimension-two in the Universal 
Mandelbrot space. Examples will be also used in other sections, where we derive (approximately) the analytical 
shape of the domains. 

3.1 Elementary domains of order p = 1 for special Mandelbrot sets 

Let us consider a Mandelbrot set for a one-parametric family 

f{x,c)=P d (x)+c (3.1) 

with polynomial Pd(x) of degree d (we do not require it to be homogeneous x d at this moment). Additive dependence 
on c-parameter considerably simplifies consideration of such families. 
For p — 1 equations (|2.1|) state simply that 

c = x-P d (x) [6 - Z) 

and for every particular choice of Pd(x) the function c(9) can be easily plotted with the help of MAPLE or 
Mathematica. Moreover, there are two important examples when even analytical solution is immediately available. 

The first case is homogeneous Pd(x) = x d , associated with the standard Zd_i-symmetric Mandelbrot sets 
Aid- The second case is generic cubic polynomial Ps(x) — x 3 + ax 2 + bx: associated family of Mandelbrot sets 
Ais(a,b) interpolates between M.2 = ^3(00, 0) and M3 — ^3(0, 0). For such interpolation one can also use a 
one-dimensional and "better" parameterized family ^2,3(0) with Ps(x) = ax 3 + (1 — a)x 2 (then M.^ = ^2,3(0) 
and M 3 = -M 2 , 3 (l))l 

3.2 Analytically solvable examples for p = 1 
3.2.1 Homogeneous Pd(x) 

For homogeneous Pd(x) = x d eq. ()3.2p converts directly into l|1.5|l : 

dx 4 ' 1 = e 1 ^- 1 ^, i.e. x = r{d)e^ 

, , (3.3) 

C = X ( 1 — X 



x(l - x d ~A = r{d)e^{l - \ e^ 1 ^) 



where r(d) = cT 3 ^ 1 and <f> — -j^. It is obvious that in this case 9 changes from to 2n(d — 1) and <fi is the right 
angle-parameter. 

Now we can use another solvable example with P3 (x) in order to deform these ideal cardioids and see how their 
order (number of cusps) can actually change in interpolation between x 2 and x 3 , see s|31 

3.2.2 Cubic polynomial 

For Ps(x) = ax 3 + bx 2 the first equation in (|3 . 2[) . 3ax 2 + 2bx = e l9 , is quadratic in x and has explicit analytic 
solution: 

-b± y/b 2 + 3ae* e . * 

x = (3.4) 

6a 

(only the "+" branch has a finite limit at a —►()). Substituting this into the second equation (|3.2p , c\ = x — P^{x) = 
—ax 3 — bx 2 + x, we obtain the analytical expression for the boundary of the root domain of the central cluster for 



2 It deserves saying that these "families" of Mandelbrot sets are somewhat artificial entities. Ai{a, b, c} is actually a 3c-dimensional 
section of the Universal Mandelbrot set, and lc-dimensional Mandelbrot sets with coordinate c are obtained if a and b are artificially 
considered as "external" parameters. Of course, one can instead take a for coordinate and b, c for parameters: no distinguished choice 
exists and all such sets should be studied on equal footing. It is nothing but a historical casus that particular sets Aid are more 
popularized than the others (and even for these particular sets the period- doubling is better known than tripling etc - despite it is in 
no way distinguished). Worse than that: the standard presentations like [5] and even the software like our favorite Fractal Explorer [4] 
implicitly exploit specific properties of these maps and produce errors in application to generic families, say, when c-dependence is not 
additive like in 113.11 1 and even if in 113.11 1 is non- homogeneous, see also introductory remarks to s|4] below. 
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arbitrary values of complex parameters a and b: 



(b T Vb 2 + 3ae ie ) (5b 3 - 9a + 3ae w =p 56V^T3ae? 



Cl = " " (3 - 5) 

Clearly, the phase transition line, separating the two regimes — 6 = 2<p (near 6 = 0) and 8 — <f) (near a = 0), 
- is |6| 2 = 3|a |. If 6 = 1 — a , see s|H it crosses the real-a line at af r = 5± ^" , i.e. a~ = 0.208712 . . . and 
a+ = 5 - a- = 4.791288 ... 

3.3 Solvable examples for p = 2 

3.3.1 Equations in case of separated c-dependence 

For p — 2 equations (|2.1| can be rewritten as follows: 

f(x) = z; 

f(z) = x; (3.6) 
f(z)f(x) = e ie 

and when 

f(x,c)=P(x) + c (3.7) 



with c- independent P(x), as 



Then c(9) can be defined from 



P(z) + z = P(x) + x; 
P'(z)P'(x) = e w 



(3.8) 



z - P(x) = x - P(z). (3.9) 



Since we did not factor out Fx from F 2 = G2-F1, these equations describe not only the (2) and (2, 1) domains, but 
also the (1) ones. The (1) domains satisfy the system (|3.8p with first equation substituted by x — z, while for the 
(2) and (2, 1) domains it should be substituted by P ^~J^ X ^ = — 1. 

3.3.2 MAPLE-generated solution for homogeneous Pd(x) 
For Pd(x) = x d the second equation in (|3 . 8|) can be solved explicitly: 

xz = dT^e iv = £ (3.10) 
where we substituted 9 = (d — 1)0. Then the first equation (|3.8p turns into 



(3.11) 



One solution, x = £/x, i.e. simply x = z = d d - 1 e 1 ^ with = % changing from to 27r, provides 



ci(0) = z - x d = x - z d = dT&e** (l - ^ e*A (3.12) 

which is our familiar eq. (|3.3|1 for the central root domain (1), with examples shown in FigO 

Remaining solutions, describing the root (2) and descendant (2, 1) domains, can be solved by MAPLE or 
Mathematica, see FigJTl In these solutions ip = <j>. No root (2) domains occur for homogeneous Pd(x) — x d , but 
this is a peculiarity of both homogeneity and p — 2: root domains (p) exist for all p ^ 2 even if Pd(x) — x d , and 
(2) domains are normally present for generic non- homogeneous Pd(x), see s|4]for examples. 
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Fi gure 7'. The first two domains (1) and (2, 1) of the central cluster, obtained from solving eq. l|3.11fl with the help of MAPLE. In 
the case of d = 2 and 3 analytical solutions are also available: see eqs. (|3,16|l and l|3.18ll respectively (and they are explicitly used by 
MAPLE). A: d = 2, f(x, c) = x 2 + c, B: d = 3, f(x, c) = x 3 + c, C: d = 3, f(x, c) = x 4 + c, D: d = 4, f{x, c) = x 4 + c. Mandelbrot set 
for f(x, c) = x d + c has Z d ~ x symmetry under rotations around the point c = 0. 



3.3.3 Analytical solutions for homogeneous P2{x) and Ps(a;) 



For d = 2 and d = 3 analytical solutions are also available. Indeed, then eas. (|3.1ip . after exclusion of solutions 
x = z. turn into 



d = 2 



X + z + 1 = x - 



1 = 0, i = -e 



i(f> 



and 



d = 3 



x 2 + xz + z 2 = x 2 + £ + 



-1, f=-e 



respectively, which are explicitly solvable quadratic and biquadratic equations. 
Then it follows that for d = 2 

1 1 / ~[ 

ar=--±-Vl-e**, 

1 1 / 



and 



°2,x{4>) = z - x 2 = x - z 2 = -1 + -e # 



is an ideal circle of radius r2,i = j centered at C2.1 = — 1, see FigJTJA. 
Similarly for d = 3 



1 1 / 2 ., 1 ., 

t = ± a / —1 e l( P ±4/1-1 e l v e 2l v 

V2 V 3 V 3 3 ' 



z = ±^=\ —1 - -e i( t> ± \ 1 + -e^ - -e 2i * 
V2 V 3 V 3 3 

and, see FigJTjB, 

C2,l{4>) = Z — X 3 = X — Z 3 

In both examples <f> = <p changes in between and 2n. 



(3.13) 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



(3.18) 
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3.3.4 MAPLE-generated solution for arbitrary cubic polynomial 

For arbitrary cubic Ps(x) the second equation in ()3.8|) is quadratic in z and can be solved explicitly. After 
substitution of this z(x) the first equation can be solved with the help of MAPLE or Mathematica and ()3.9j) 
produces the final answer. This is how Figs are obtained. 

4 The first two elementary domains in interpolation between M.2 and 

After equations are solved, we can turn to description of solutions. 

For particular homogeneous polynomials Pd(x) — x d we obtained the well known shapes of central domains in 
Md Mandelbrot sets, see FigUJ - only this time this is not a result of computer simulations by Fractal Explorer, 
the shapes are now obtained as solutions (sometime even analytical) to algebraic equations (|2.1|) . 

Even more interesting is the possibility to study quantitatively interpolation between the M.i and M3 sets. 
So far only qualitative description was known [2], and the usual computer simulations fail. Such simulations [5] 
are often based on the study of the sequence / on (c), i.e. the orbit of x cr = 0. Interior of the Mandelbrot space, 
i.e. the black regions in FigsHHU is assumed to consist of all functions / where this sequence is bounded and does 
not go away to infinity. However, this assumption is not always true and then this simple algorithm fails - and 
Figgis the first example. The reasons for failure can be different: from x cr ^ to attraction of unstable orbits to 
finite, rather than to infinitely remote ones. There is a strong need to cure this problem and make a modification 
of Fractal Explorer which would treat properly any kind of Mandelbrot setd 

4.1 Particular Mandelbrot sets ^2,3(0.) for the families ax 3 + (1 — a)x 2 + c at different 
values of a 

At a = the Mandelbrot Set acquires its standard form, Fig[TJ and its first two domains, (1) and attached (2, 1), 
are shown in Figs[TjC and [JJ A. 

4.1.1 Vicinity of the Mandelbrot Set: small \a\ 

However, as soon as a infinitcsimally deviates from a = 0, it gets clear that Fig[T]has a twin: an exact copy of the 
same shape and size, but with opposite orientation - a mirror twin, - located at infinity of the complex-c plane. 
As \a\ grows, the twin moves closer, and Figs [8] and [9] show its location at a = ±1/10 (the sizes of the domains 
are practically the same as in FigQ]- just the scale of the picture is different, because the twin of FigQ] is still 
far away). Moreover, it appears that additional mirror pair of (2)± domains - roots of two more clusters - were 
hidden at infinity of c plane and are now located in between the two root domains (1)± for positive a > and on 
the opposite sides of those for negative a < 0. Since for small a these domains are tiny as compared to the (1) and 
(2, 1), they can be easily overlooked, therefore one of them is marked by a circle and shown in a bigger scale in a 
separate picture at the right low corner. Clearly this root (2)_ domain has cardioid shape and is exact copy of the 
root domain (1)-, only smaller. In fact it has a (4, 2)_ domain attached to it in exactly the same manner as (2, 1)_ 
is attached to (1)- - it is not shown, because we explicitly construct only domains of orders v = 1 and 2. For 
interested reader we add also slices of the Julia sheaf: show behavior of the orbits in the x spacqj with the change 
of c, which becomes more and more interesting as we go far from the "pure" points a — and a = 1. The problem 
is that Julia sheaf is embedded into a 4#g? = 2cd space, with complex x and c, and can not be shown in full, even if 
a is fixed. Therefore different sections and projections are presented, 2^- and 3fl-dimensional. 3i?-dimensional are 
especially informative, but only when presented on computer screen, where they can be rotated and regarded from 
different angles. This advantage is lost in the printed version of the text, but one can either use simple MAPLE 
programs, collected in s[T2] below or directly look at the results in [BJ. 



3 In the absence of such modification we had to make use of various pictures, which are at best qualitatively, but not fully correct: 
this is the case with Fig|4]in this paper and with numerous Figures in [2], including even the picture at the cover of that book. Below 
in this section we provide much better views of the 2-parametric section of the Universal Mandelbrot Set, these pictures will be fully 
correct, but instead only order- 1 and 2 domains will be shown. 

4 Since f(x) is cubic, there are three order-1 orbits and up to three branches will be seen in the pictures. Since G2(x) = F2(x) / F\{x) 
has degree 6 in x, there are 6/2 = 3 orbits of order p = 2 and up to six branches will be seen in the pictures. 
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Figure 8: The picture in the left upper corner represents the Mandelbrot set for the family ax 3 + (1 — a)x 2 + c with a = 1/10. This is 
the picture in complex c plane and shown are only domains, associated with the order p = 1 and p = 2 orbits. The orbits themselves lie 
in the complex x plane and form the Julia sheaf over Mandelbrot set. Julia sheaf itself is a led complex variety, embedded into a 2(jd 
complex space and can not be shown in an ordinary drawing. Instead two different views, one 2^-dimensional, another 3fl-dimensional 
(at fixed Re(c) = C(a)) are shown in the low left and upper right corners respectively. As all other three dimensional sections in this 
paper, it can be rotated and looked from different angles: this clarifies the pattern a lot, but can be done only on computer screen, see 
ref. [6], (In particular, there are NO intersecting orbits in this section - the seeming intersection is an artefact of the drawing, resolved 
by rotation of the picture). Dilute lines represent the three order-1 orbits, while dense lines — the six order-2 orbits. In 2nd picture 
solid segments show stable orbits: those of order 1 are stable inside the (1)± domains of the Mandelbrot set, those of order 2 - inside 
the (2, l)± and (2)± domains (the last two are two small to have the corresponding solid segments seen in our pictures). Enlarged 
picture of the (2)_ domain - a root of a new cluster — is shown in the low right corner, and it is clear that it is an exact diminished copy 
of the root (1)- domain. (2)+ is its exact mirror copy, in accordance with Z2 symmetry of the Mandelbrot set w.r.t. the vertical line 
Re(c) = C(0.1) = —8.4. Shown are also the roots c = d\ } 2 of discriminants D\ and D2 and c = ri2, c = T24 of the resultants R12 and 
R24, (they lye at intersections of vertical lines with the real-c axis). According to [2], the last two are the crossing points of the orbits 
of orders 1 and 2 and 2 and 4, define merging points between the domains (2, 1) and (1) and between (4, 2) and (2, 1) respectively, 
and thus define the stability segments of orbits of orders 1 and 2. Similarly, discriminant zeroes are intersection points of the orbits 
of the same order: p = 1 with p = 1 and p = 2 with p = 2. 2nd view in the low left corner is in fact a section of the c plane with 
given Im(c) = and projection on the Re(x) plane. Accordingly, when two real-valued orbits intersect and become complex-valued, 
they remain shown in the picture, but since they are complex conjugate, two lines are projected onto one - this should be taken into 
account in analysis of the figure. 
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Figure 9: Analogous picture for negative real a = —1/10. The only essential difference from Fig|8]is that the two root domains (2)± 
are not in between the two (1)±, but on the opposite sides of those. In other words, at a = the (2) domains pass through (1)-, so 
that (2)-|_ re-appears from c = +oo, while (2)_ returns back from c = — oo, but slower than (1)—, see also Fig llOl Because of this, 
when \a\ increases further in the direction of negative a, the overlap will occur between the domains (2, 1)_ and (2, 1) + , unlike in the 
positive-a case, where (2)_ and (2)+ will be the first to meet. 
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Figure 10: Behavior in the vicinity of a = of the root domains (1)— and (2)±, denoted respectively by circles and solid lines???. As 
a — > 0, these three domains, together with their entire clusters, travel to infinity in the complex-c plane, therefore the picture is drawn 
in coordinate 1/c (only real values of c are plotted). The Mandelbrot Set from Fig[T] including the central root domain (1) + , stays in 
the vicinity of c = and is not shown in this picture. Behavior of the (2)_|_ and (2)_ domains is somewhat different: the former one 
re-appears from infinity at the opposite end of the Re(c) axis, while the latter one returns from the same end, only exchanges positions 
with the (1)— domain — in full accordance with Figs[8]and[9] All the domains and in fact the entire clusters shrink to a single point at 

a = 0, the reason for this is the choice of a very singular map in homogeneous coordinates: (^x, y\ > (ax^ + bx 2 y + cj/ 3 , y 3 ). Such 

singular behavior at a = would be smoothed and become similar to bifurcations at finite a, shown in forthcoming pictures, if y s is 
substituted by generic cubic polynomial. See [3] for related considerations. 
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All Mandelbrot sets ^2,3(0) possess a discrete symmetry under reflection w.r.t. the vertical line 



3a \ 9a 

with b = 1 — a, which is lifted to entire Julia sheaf over .^2,3(0): 

c = c—C(a) — ► — c 



Re(c)=C(a)=-^-(l + ^-\ (4.1) 



(4.2) 



For example, the equation Fi(x,c) = for the first-order periodic orbits is obviously symmetric, since Fi(x, c) = 

ax 3 + bx 2 + c — x = ax 3 — ^1 + |^ x + c. In accordance with this symmetry, the (1)_ domain - the twin of the 

(1)+ domain, centered at c — 0,- is the mirror-reflected cardioid with center at c = 2C(a). Similarly, the two next 
root domains (2)± are centered at c 2 + = — - + 2 + 0(a) and c 2 _ = 2C(a) — c 2+ , see s.6.4.5 of ref.[5]- Since c 2+ is 
odd function of a for small a, while Ci_ is even, it is clear that they exchange order when a goes through zero - in 
accordance with what is shown in Figs 151 151 and ITD1 



4.1.2 Overlapping domains 

As |a| increases, the two domains (1)± move closer. The speed of approaching is somewhat different for positive 
and negative a. At some stage of this movement two different clusters unavoidably meet. There are, however, two 
sorts of meeting: overlap and collision. For above-explained reasons it is still difficult to analyze the behavior of 
entire clusters. Approximate results about clusters (or, better, possible approach to their future derivation) will be 
discussed in the last sections [STTTOl of this paper. Exact results to be considered right now concern meetings of the 
low (1 and 2) order domains. Overlap of these domains takes place soon after it happens to the upper (p = 00) 
leaves of the clusters, while collision can start at the p = 00 leaves, but can also begin at the low-p level. 

Overlaps and collisions of particular domains occur when zeroes of the corresponding resultants collide, i.e. are 
controlled by zeroes of double-resultants, like c-discriminants of the x-resultants, listed in the following table (italic 
lines are quotations from MAPLE program, boldfaced are real- valued roots, belonging to the segment < a < 1). 



Two zeroes of Dl merge at a = — 1±l £^ . 

discrim(discrim(Fl, x), c) = 96 * a 2 -I- 112 * a 3 + 96 * a 4 + 48 * a + 48 * a 5 + 16 + 16 * a 6 
Two zeroes of R12 merge at a = 5± ^i and at a = — 2 ± \/3. 

Only one of these points, a = 5 ~^~ — 0.20871215 . . . belongs to the interval < a < 1 on a real-a line. 
discrim(resultant(G2, Fl, x), c) = 16 * a 10 * (a 2 + 4 * a + l) 2 * (a 2 — 5 * a + 1) 
Two zeroes of D2 merge at a — 4 ± \/15. 

discrim(^J discrim(G2 1 x)/ R12 j ~ (a 2 — 8 * a + l) 3 

Two zeroes of R24 merge at 

a = 0.1487496031 ± 0.03597329725?, a = 6.351250397 ± 1.535973297*, 

a= -4.485250968,-0.2229529645, 0.1163898556, 8.591814077= 1/0.1163898556, 

a = 0.1497297977, 6.678697327 = 1/0.1497297977, 

a = 0.5857864376 ± 0.8104654524i, 

_ _ 9±y / 65 
a — 4 . 

The four series correspond to zeroes of the four factors in 

discrim x), cj (2 * a 4 — 26 * a 3 + 93 * a 2 — 26 * a + 2)- 

•(a 4 - 4*a 3 -39*a 2 -4*a + l) 2 (a 4 - 8*a 3 + 10*a 2 -8*a + 1) 2 (2 * a 2 - 9 * a + 2) 4 



Derivative dd^[ /da for the zero d^ of discriminant D±, which defines the position of the cusp in the (1)- domain, 
defines the speed of motion of the cluster (1)- with the change of a. Since for small a all the clusters are diminished 
copies of the central one in Fig[TJ with all the same proportions one can actually predict what happens to clusters 
from the data about their root domains. 

The first event to happen on our way from a = is overlap. If a > this is the overlap of domains (2)+ and 
(2)_, while if a < it is that of (2, 1) + and (2, 1)_. The two domains of order 2 overlap when the two zeroes of 
i?24 coincide. If we move from a = along the real-a line this first happens at a = 0.1163898556 ... if a > and 
at a = -0.2229529645 . . . if a < 00 FigsIHl and [T2l show Mandelbrot sets soon after these points are passed. 

5 In fact, as explained in the previous paragraph, we can approximately find the values of a c i ugter , when the overlap of the clusters 
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Figure 11: Domains of the orders p = 1,2 of the Mandelbrot set -M2,3(a) (domains of orders p = 1,2) at a = 0.117, immediately 
after the two root domains (2)-|- touched at a = 0.1163898556 . . .. Now they overlap: for c in close vicinity of c = —6.24 there are two 
stable order-2 orbits at once. However, as clear from the Julia-sheaf pictures around, nothing special happens to the orbits themselves. 



It is clear from these pictures, that when overlap occurs, nothing interesting happens to the orbits - and this 
is what makes overlap different form collision, when intersection of orbits takes place, see below. Overlap means 
simply that two (or more) different orbits are simultaneously stable at the same value of c: in this case these are 
two order-2 orbits. When overlap increases, it involves the (1)± domains and at the same values of c can coexist 
two stable orbits of other orders: 2 and 1 fFig fl"3]) and 1 and 1 (FigfM]). For a < — 1 the further increase of \a\ leads 
to diminishing of the overlap: the story repeats in the opposite order, the overlap picture for a = —4 ressembles 
that for a = -0.25 (shown in Fig fig]), for a = -4.(4) - that for a = -0.225 (shown in Fig|T2"|). and after that the 
overlap disappears. 

The reason for reversed evolution is that the family M2, 3(d) and even its Julia sheaf has a discrete "symmetry" 
w.r.t. inversion of parameter a — > 1/a: 

x(a~ ,6) = —ax(a,9), 
c(a-\6) = -ac(a,6) (4.3) 

This symmetry complements the (|4.2[) of particular Mandelbrot sets at fixed a, and it allows to consider only the 
variation of a within interval \a\ < 1, all Mandelbrot sets outside this interval are exact rescaled copies of those 
inside. 

More interesting things are taking place for positive a. 



occurs. From Fig[T] we know, that the total size of the cluster is 1.65 times bigger than the size of the root domain, and the 
latter size is nothing but the difference 0I2 — r24, we get a rough estimate: a^ — a c i U3ter PS 0.65(a co — a^ a ) where and a co are 
moments when ?"^(a [ j ) = r^ia-do) an d rf^(a C o) = ii^(a co ), i.e. a^o = 0.1163898556 . . . and a co = 4 — %/l5 = 0.12701665 . . .. Then 
a cluster ~ 1.65a,j — 0.65a co PS 0.109. It is assumed that the speed of motion of the cluster with the change of a and the cluster's size 
are approximately constant, appropriate corrections can be easily taken into account. 
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Figure 12: Domains of the orders p = 1, 2 of the Mandelbrot set A^2,3(a) at a = —1/4 + 1/40 = —0.225, soon after the two descendant 
domains (2, 1)± touched at a = —0.2229529645 . . .. Two order-2 orbits are simultaneously stable in the close vicinity of c = —0.88. 
Again, nothing special happens to the orbits themselves when the overlap occurs. 
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Figure 13: Domains of the orders p = 1, 2 of the Mandelbrot set A^2,3(i) at a = —1/4. The overlap increases, and now the descendant 
(2, l)± domains overlap not only each other, but also the twin parent domains. This means that in the vicinity of c = —0.65 

there are two coexisting stable orbits of order 2, while in the vicinities of c = —0.5 and c = —0.8 coexisting are stable orbits of orders 
2 and 1. Nearly cusp-like behavior of the orbits in Julia-sheaf view in the upper right corner occurs in the vicinity of the zeroes ri2 
of the resultant R12, where the orbits of order two cross those of order one and the (2, l)± domains are attached to their parent (1)± 
domains. This singular behavior has nothing to do with the overlap of the domains. 
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Fi glire 14: Domains of the orders p = 1, 2 of the Mandelbrot set A / i2,3(^) at a = — 1. The overlap increased even further as compared 
to Figs ll2l and 1131 Descendant domains (2, 1)± now passed through the root ones (1)=f and are no longer involved in the overlap. 
Instead overlapping are the root domains (1)± and two stable orbits of order 1 coexist in the vicinity of c = 0. If (1)± domains 
continued to move in the same direction with further increase of \a\, the two order 1 orbits would cross (when the two zeroes df of 
discriminant D\ merge), and overlap would finally end in a collision - like it happens for positive values of a, see Fig |16l However, the 
merging of D\ zeroes occurs at complex values a = — ^\ > ant ^ at a = — 1 the overlap is the biggest possible for real negative values 
of a. As a decreases further (and \a\ grows) along the real-a line, the (1)± domains start to move in the opposite direction, see Fig |15l 
We pass again through the overlap patterns like Fig |13l (at a = —4) and Fig |12l (at a = —4.(4)) and finally come back to the no-overlap 
pattern like Fig|9](at a = —10) - all this in accordance with the symmetry 1 14. 31 . 
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o=-2 




Fi glire 15: Domains of the orders p = 1, 2 of the Mandelbrot set A^2,3(^) at a = — 2. The overlapping clusters start to diverge after 
maximal approach at a = — 1, stopping short from collision of two zeroes of discriminant D\, which would unify two clusters into a 
single one This actually happens at complex values of a = 1 ^ v ^ , and we will enconter this unified cluster in our travel along real-a 
axis, but at positive values of a. 
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4.1.3 Colliding domains 



We left evolution in the positive-a direction at the stage of FigfTT| when overlap of the two root domains (2)± just 
occured. In variance with the case of negative a, this time as the overlap increases the two zeroes di of discriminant 
Z?2, defining positions of the cusps of these two domains, coincide at a — 4 — y/l5 — 0.12701665. . . and a new 
phenomenon takes place. The two stable order-2 orbits (they were simultaneously stable in the overlap region) 
cross each other, and the pattern of orbits around crossing is pretty sophisticated, see Fig[16] Most interesting is 
preservation of two small overlap sections, where two different stable order-2 orbits continue to coexist, but exhibit 
non-trivial monodromy under a travel in the complex-c plane around the cusps at zeroes of discriminant Di- 



4.1.4 Colliding clusters 

After collision of two root (2)± domains and formation of a unified cluster with the root (2), the two other clusters, 
growing from the (1)± roots, continue to move towards each other and soon collide with the (2) cluster, sandwiched 
in between them. Now this is indeed a collision, not just overlap, and, in variance with collision of the (2)± domains 
it now originates at the highest leaves (at p = oo) rather than at the root domains. Full description of this process 
is impossible with the knowledge about the p = 1, 2 orbits only, thus our illustrations will be necessarily incomplete. 
Still, a lot is seen even with these limited tools. 

Immediately after collision of the (2)± domains in Fig ll6l they begin growing and soon become comparable in 
size with the (2, 1)± domains, belonging to approaching clusters. Even earlier the overlapping region inside (2) 
shrinks down and disappears. Finally, when the zeroes ri± of the resultant R24, marking the closest points of the 
(2) and (2, 1)± domains, coincide, collision wave, going down from the upper leaves of the clusters, reach the p = 2 
level: cluster collision gets seen at the level of our consideration. This happens at a = 0.1497297977 . . ., see Fig l 171^1 
FigsfTBl and [T9l show Mandelbrot sets soon after that and a little later, when continuing approach of (1)± domains 
(which are now two roots of a single cluster!) starts pushing the unified (2, 1) domain outside of the region between 
them. This push-away process leads to the next bifurcation at a = 5 ~^"" = 0.20871215 . . . where the two zeroes 
r\2 of another resultant R12 coincide, marking collision of the (1)± domains: collision wave reached the p = 1 level. 
At this moment the (2) cluster is ripped into two disconnected pieces, see Figj20] Clearly, just the same push-away 
and ripping processes took place with all the higher-order domains (2 fe , 2 fc ~ 1 , . . . , 2, 1) in between the moment of 
clusters collision till it reached the p = 2 level at a = 0.2087 . . . 



4.1.5 Mandelbrot sets with the topology of Ai^ (in the vicinity of a = 1) 

The further evolution of ^2,3(0) with increasing a consists mostly of continuous deformation of the shape of 
emerged unified root (1) domain with two cusps: from a bone-like region in Fig |20l it grows into nearly oval one 
(deviating from oval only near the cusps) in Figs l21l 1221 and finally at a — 1, when the cusps extend their region of 
influence, acquires the standard form of Fig l23l familiar from Figf2] 



Still, while nothing equally drastic happens after the cluster collision, evolution is not quite event-less. In above 
pictures one can see that the overlap regions in the (2, 1) domain(s) appear and disappear, signaling about the 
motion of orbits in Julia sheaf with the change of a. In particular, an interesting inside-out reshuffling is shown in 
FigETJ Moreover, in Fig [22] one can see that overlap occurs even between the (1) and (2, 1) domains. We emphasize 

6 Like in footnote [5] one can try to estimate the moment of clusters collision. However, this time the clusters shape deviates 
considerably from that in Fig[T] therefore such estimate is less reliable. 
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Figure 16: Domains of the orders p = 1,2 of the Mandelbrot set at a = 1/8 + 1/224 = 0.12946(428571). The two root domains 
(2)_)_ and (2)_ just collided at a = 4 — \/TE = 0.12701665 . . . (after a period of overlap, originated at a = 0.1163898556 . . ., see FigfTTJ, 
and their clusters are already united into a single new cluster, so that the union of (2) = (2)_|_ U (2)_, (2)+ n (2)_ 7^, is now the new 
root domain. Arrows in the picture at the low left corner show the action of evolution x — > f(x) on the points of the order-2 orbit. 
Enlarged picture in the left upper corner shows in more detail the new root cluster (2). It has two narrow regions of self-overlap along 
the Re(c) = Re(d2) vetical axis, where two different order (2) orbits are simultaneously stable. In all other points of the (2) cluster 
only one order 2 orbit is stable, but non-trivial monodromy occurs when we go around the cusps at = Ke(d2) i Im(d2) (points D): 
if we pick up the single stable order-2 orbit, say, at c = Re(d2) (point O) and carry it into the upper overlap region from the right 
(counter-clock-wise), we obtain one of the two orbits, stable in that region, but if we carry the same orbit into the same region from the 
left (clock- wise) , we obtain another stable orbit. If we continue to carry our orbits in the same direction and leave the overlap domain, 
the orbits loose stability. In other words, the two stable orbits in the overlap domain are permuted when carried around the cusp, and 
each orbit continues to be stable outside the overlap domain only if it is carried away in one out of two possible directions: either to 
the right or to the left. Another end of the overlap domain (point E), which is on the boundary of the (2) cluster is not a cusp and not 
a singularity. As any point on the boundary of an elementary domain it has in its infinitesimal vicinity an infinite number of zeroes 
of various resultants R2.2m, where the stable order-2 orbit exchanges stability with some order 2m orbit, see [2]. However, the set of 
relevant {m} changes irregularly as we move along the boundary, and accordingly, for the point E this set depends irregularly on a. 
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Figure 17: Domains of the orders p = 1,2 of the Mandelbrot set Ai2,:i(a) at a = 3/20 = 0.15. The (2) domain just merged with 
the (2, l)± at a = 0.1497297977 ... to form a single descendant (2, 1) domain. This unified domain has a characteristic four-sausage 
structure, remembering about its recent formation from four distinct elementary domains (2)± and (2, 1)±. 
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Figure 18: Domains of the orders p = 1,2 of the Mandelbrot set M 2 ,3( a ) at a = 1/7+1/84 = 0.15(476190). The shape of the unified 
(2, 1) domain evolved from the four-sausage to a two-sausage one: the memory about the difference between the (2) and (2, 1) domains 
is almost erased, but distinction between + and — is still preserved. 
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Figure 19: Domains of the orders p = 1,2 of the Mandelbrot set M2,3(a) at a = 1/5 = 0.2. The (2,1) domain is being pressed 
away by approaching (1)± domains, which are now the two roots of a single cluster. Note re-appearance of the narrow overlap regions 
inside the (2, 1) domain and non-trivial monodromy of the order 2 orbits when they are carried along a circle in the Mandelbrot set, 
surrounding one of the cusps <i*. 
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0=1/5 H/80 




Figure 20: Domains of the orders p = 1, 2 of the Mandelbrot set M 2 ,3<a) at a = 1/5 + 1/80 = 0.2125. The two root domains (1)± 
just collided at a = ''~?f^ = 0.20871215 . . . and formed a single root domain (1). It has a typical bone-like shape, with two pieces of 
the former descendant (2, 1) domain attached at the merging point. Overlap regions inside these (2, 1)± fragments are now pretty large 
and cusps are not well seen, still the picture in the low right corner shows that they are still there. 3rcI view in the low left corner is 
taken at Im(c) = 0, where only order p = 1 orbit can be stable inside the root (1) domain. 
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Fi glire 21; Domains of the orders p = 1,2 of the Mandelbrot set A^2,3(i) at a = 1/4. In this picture we are near the point 
a = = 0.23443556 . . . where inside-out reshuffling of the (2, 1) domains takes place (two zeroes r24 of R24 coincide at this 

transition point). Overlap regions, seen in Fig l20l grew up to the full size of the (2, 1) domains, and boundaries of the regions changes 
roles with boundaries of the domains. Later the former boundary domains will get closer and form new small overlap regions in Fig |22l 
that will later disappear at a = 0.42 . . . 
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Figure 22: Domains of the orders p = 1,2 of the Mandelbrot set M2,z(a) at a = 1/3. The root (1) domain acquired an almost 
ellipsoidal form (outside the cusp regions). The interesting phenomenon, seen already in Fig |21l is appearance of additional overlap 
regions, marked by arrows in the upper left corner: between the root (1) and descendant (2, 1)± domains, where order-1 and order-2 
orbits are simultaneously stable. This time there is no interesting monodromy: when we leave the overlap region in the direction of (1) 
the order-1 orbit remains stable, when we go deep into (2, 1)±, the stable one is an order-2 orbit. 
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Fi gure 23 ! Domains of the orders p = 1, 2 of the Mandelbrot set A42,3 (a) at a = 1. This is the central part of the standard Mandelbrot 
set M3, shown in Fig[2] All overlap regions disappeared, the (1) domain is ideal cubic cardioid, other elementary domains acquire 
nearly cardioid shapes. 
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once again, that no bifurcations (phase transitions, orbit crossing or reshuffling) are associated with the overlaps, 
still they affect the shape and even the very presentation of the Mandelbrot set (when overlaps exist, it is not 
a clever idea to draw it all in black, like we did in FigsQ][4|) and in fact this is a signal that the phase portrait 
gets richer: a non-trivial pattern of attractors and repulsers occurs, nothing to say that the vicinities of unstable 
orbits are not necessarily attracted to infinity, as implicitly assumed in some algorithms, mentioned in the first 
paragraphs of s[4j 

Events, encountered in the evolution of Mandelbrot set ^2,3(1) from a = to a = 1, i.e. in interpolation 
between FigsQ]andEl are collected in the following table: 
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a 


typical feature 


picture 


a < 






a = 


the standard Mandelbrot Set M2 


FigQ] 


< a < 0.1164. . . 


two root domains of type (1); 
two descendant domains (2, 1), attached to them; 
two isolated root domains (2); 
two (4, 2, 1) domains and two (4, 2) domains, 
attached to (2, 1) and (2) at real values of c 


FigE 


a = 0.1164. . . 


projections of two domains (2) meet, 
i.e. two zeroes of R24 coincide, 
responsible for stability of two different order-2 orbits 




0.1164. . . < a < 0.1270. . . 


projections of two domains (2) overlap; 
two stable order-2 orbits coexist 
in the region of overlap 


FigE] 


a = 4- v / l5 = 0.1270... 


cusps of overlapping domains (2) merge, 
i.e. two zeroes of D2 coincide 




0.127. . . < a < 0.130. . . 


overlap of the two domains (2) 
splits into two isolated components 


FigED 


a = 0.130 . . . 


overlap region shrinks down 




0.130 ...< a < 0.1497 . . . 


only one unified root domain (2) exists 




a = 0.1497. . . 


(2, 1) domains merge with the (2) domain, 
i.e. two pairs of zeroes of R24 coincide, 
each pair responsible for stability 
of the same order-2 orbit 


Figna 


0.1497. . . < a < 6.6800. . . 


only one (2, 1) domain exists; 
no (2) domains 


Figffl 


a = 0.1875 . . . 


coexisting stable order-2 orbits re-emerge 


Figfl9l 


o= 5 ~^" = 0.2087... 


two (1) domains meet and 
the (2, 1) domain splits into two, 
two zeroes of R12 coincide 


Fig®] 


0.2087.. . < a < 0.42. . . 
0.2087 . . . < a < 4.7916 . . . 


overlap regions, where two stable order-2 orbits 
or order-2 and order-1 orbits can coexist; 
one (1) domain; 
two attached (2, 1) domains; 
no (2) domains 


Figsl2T1fel22l 


a = 0.42. . . 


overlap region shrinks down 




a = 1 


the standard Mandelbrot set .M3 


I*igsl23lfel2l 


a > 1 








1* III 



Figure 24: Domains of the orders p = 1,2 of the Mandelbrot space M2,3 (a) at a = 5. This picture is direct analogue of Fig |19l 
and serves as an illustration of the symmetry property 14.31 1: that the same sequence of bifurcations happens to the Mandelbrot set 
:t(a) on the way from a = 1 to o = -1 along the real-a line through a = oo as on direct way, presented in Figsl8 H23l 



In s !4. 1.21 we briefly discussed what happens beyond the realm of this table: for negative values of a. The 
evolution of M.2,3(a) can be also continued to the region where a > 1 (where b = a — 1 < 0). This evolution 
appears to be reverse of what we already considered: the Mandelbrot set of FigfSJat a = 1 passes through the same 
stages of FigsJ22 | [21] , [201 (at a = 3, 4, 4.706 . . . respectively) and so on. In particular, at a = 5+ ^ T = 4.7912878 . . . 
the single root (1) domain is split into two, while two descendant (2, 1) domains merge into one. For illustration 
we show in FigJM|the counterpart of FigJTni The full picture will be shown in s !4. 21 below. 



4.2 First steps towards UMS 

Let us now return to interpretation of Mandelbrot sets as sections of a single Universal Mandelbrot Set (UMS). 
It implies that all what we observe about particular collection of Mandelbrot sets, like our one-parameter family 
■M2, 3(0)1 can be re-interpreted as result of particular view on one and the same solid structure: variation of patterns 
is result of the changing view, the structure is always the same0. 

For example, the entire collection of the p = 1 domains in Figsl51l2~4l which consist of a single component for a 
in between 5=F ^'" and split into two components outside this segment, see Fig[25lA, can be alternatively described 

7 We can not avoid stressing analogy with the well-known projection approach to integrable systems, see [9] 
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Figure 25: A. The shapes of p = 1 elementary domains for different values of a in the family ax 3 + (1 — a)x 2 + c, as given by 
eq, l|4.4j l, B. The same (qualitatively) pictures arise in the sections of a single cardioid-like cylinder by different parabolic-like sections 
[x — (z 2 /2 — s) = for various s) with a complex-c plane. C. The cylinder can be even made circlic??? with the help of the 
circle-cardioid relation shown in Fie l26l The corresponding sections are now generic conies (quadrics???), not obligatory paraboloid. 
The true sections behind ea.l l4.4l l are cubic, see 1 14. 5t . not a quadric, which are a little less convenient to draw. 



as the image of a single cardioid-like domain in the slices, evolving with the change of the section in Fig[23lB. 
Moreover, in this approach one can even start from an ordinary circle, not from a cardioid, see Fig |26l In fact these 
pictures are nothing but approximate drawings, attempting to capture the properties of exact formula 



_b_( 26*\ / 2b 2 e t9 \ Vb 3 + 3ae' e 

~3^ V + 9a / ± V + 9«" ~~ ~3~/ 3a ' ^ ' ' 

just now we interpret it is an evolution of (degenerate) elliptic mappings 

(c — co) 2 = k(u — u\)(u — U2~)(u — Us), with a — dependent parameters Co, k and Ui] actually U2 = U3 (4.5) 

of a complex-M plane with a unit circle on it into a complex-c plane, where the image of the unit circle looks 
differently: like our p = 1 domains, evolving and even bifurcating (splitting and merging) under the change of the 
mapping. 



At this early stage of investigation of UMS it is unclear, what is its best and simplest possible representation. 
In particular, nothing as simple as Fig ]25l is immediately available when the order-2 orbits are taken into account. 
Therefore, instead of playing with different realizations, we show in FigsETl [28] the most straightforward views of 
the Universal Mandelbrot Set, directly in the a — c coordinates. Unfortunately, only 3i?-dimensional section of the 
full 2c-dimensional pattern can be drawn in a picture and it can be rotated (what is very informative!) only on 
computer screen, see [5] for details about this option. 
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Figure 26: Cardioid cylinder of the previous Figure l25l can be substituted by an ordinary rotation-symmetric cylinder at expense 
of a more sophisticated slicing. Cardioid as a square of a circle. Analytically this correspondence is represented by c = t| (l — or 
1 - 4c = (u- l) 2 . 

In Fig[29]we give some more examples of 2c-sections of the UMS. In particular, we demonstrate that topologies 
of this section can be very different and they can be investigated by already available tools. As a small illustration, 
the chain of pictures in FigfJUB shows how a loop in particular section of the UMS can be (on Figf^HA) contracted. 
Of course, we are very far from calculation of homologies of the UMS, but the way is already open. 



5 Small-size approximation (SSA) 

We now switch from transparent exact results to subtle approximate methods. 

5.1 On status of the SSA 

The shape of elementary domains can be also considered in the "small-size" approximation (see s.4.9.3 of [5] and its 
less accurate predecessor in 8 and many other text-books). In SSA we expand all functions of x and c in powers 
of their deviations from the critical (or, simply, mean) values and leave only the first three (constant, linear and 
the next) terms in these expansions. In what follows x cr = 0, but c cr will have different values. "Next" normally 
means quadratic, but for homogeneous polynomials Pd{x) = x d (giving rise to Zd_i-symmetric Mandelbrot sets) 
it will actually mean x d . 

SSA would be very natural, if typical deviations of x and c from the mean values were small. However, while 
it can seem reasonable for the study of elementary domains - except for the first few, they are indeed pretty 
small in the c-plane (most are not even seen in Figs[TH3]),- actually this assumption is wrong: as we already know 
from Figs[S][Ml the x- variables in solutions to eqs. (|2.ip are not small. At today's level of knowledge justification 
of SSA comes only a posteriori, say, from comparison with experimental data in sJSJ) . SSA seems adequate for 
phenomenological description of experimentally observed [5] self-similarity (fractal or scaling) property of the 
Mandelbrot sets, but no clear theoretical reason for this adequateness is known. The algebro-geometrical approach 
of [2] only adds to the mystery: the (experimentally) obvious scaling properties of the universal discriminantal 
variety call for clear conceptual explanations. 

In any case, today SSA is the only available approach to evaluation of Feigenbaum indices and other charac- 
teristics of elementary domains when their order p — ► oo. Surprisingly or not, it does a very good job in this field: 
see slHl below. 

5.2 SSA for the Mandelbrot Set 

The first step of the SSA in application to the Mandelbrot Set (i.e. to the family f(x;c) = x 2 + c) is to expand 

F P (x) = f p (c) -x + x 2 lp (c) + 0{x A ) (5.1) 
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Fi glire 27i Collection of order-1 and -2 domains of Mandelbrot sets A^2,3( a ) f° r the family f(x) = ax^ + (1 — a)x^ + c with various 
values of a in the region 0.10 < a < 0.28, where most bifurcations are taking place, represented as slices of a single 3r<1 entity, which 
can be also considered as a 3rcI section of UMS. Horizontal is the plane of complex c, vertical is the line of real a. This picture provides 
a concise summary of all the properties, described throughout s l4.ll The lower part of the figure represents separately the p = 1 and 
p = 2 domains: all pictures are also shown from different angles, what can help to appreciate the beauty of the structure. 



39 



If 




Re c 



Figure 28: The same collection of Mandelbrot sets A^2,3(«) for the family f(x) = ax 3 + (1 — a)x 2 + c with all real values of a 
(including those represented in Figs l8l24t . Only central domains of order 1 are shown. 
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Fi gure 29: Examples of other 3^rf sections of the Universal Mandelbrot Set. Only central domains of order 1 are shown. 
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It is used here that that f'(x = 0) = 0, otherwise expansion in powers of x should be substituted by that in powers 
of x — x cr , where f'(x cr ) = 0. 

The second step is to substitute (|5.ip into (|2.1jl . This gives: 

2x lp (c) = e l<t> (5.2) 

and 



2/ p (c) 7p (c) = 2x lp (c)(l - x lp {c)) *P e« 



e*\ (5.3) 

Eq. (|5.3p provides the answer: the l.h.s. depends on the shape of the function f(x;c), i.e. on c, and (|5.3p describes 
a curve c((j>). 

Actually this curve can have many disconnected components, and the third step is to consider a particular 
component, surrounding a particular root c p of the Mandelbrot function f p (c) — f(x cr ,c): 

f p (c p ) = (5.4) 

Expand (15. 3p around c p : 

c = c p + a. 



fp(c) = f P (c P + a) = f p a(l + Aa + 0(a 2 )j (5.5) 

7 P (c) = 7p(c p + cr) = 7 P ^1 + ^ CT + 0(er 2 ) j 

From now on we denote through f p and 7 P the values of the corresponding functions at c = c p : "f p = j P (c p ) etc. 
Substituting ()5.5|) into (|5.3[) . we get: 



with 



and 



2/ p 7p 



(5.6) 



(5.7) 



ep = -^-(4 + ^] (5.8) 
/p7p V 2 /p Tf/ 

Eq. (|5.6p is our final SSA answer for the shape of the elementary domain of the Mandelbrot Set, surrounding a 
point c p . We see that the complex- valued r p defines the size and orientation of the domain, while its shape is fully 
controlled by the value of £ p : for <C 1 we get a cardioid, while for |£ p — 1| <C 1 it turns into a circle. 



Thus our problem is reduced to: 

- the check of accuracy of the small-size approximation: we do this in s[6] by comparing the values of r p , 
predicted by (|5.7|) with their actual values for the family f(x; c) = x 2 + c, measured with the help of the Fractal 
Explorer j¥j or defined from the roots of the relevant resultants; 

- evaluation of parameter £ p with the help of (|5.8|) : in s[7]we show that indeed in the small-size approximation 
£ p = 1 for elementary domains, which are not roots of any clusters (i.e. are descendants of some lower-level 
domains); 

- demonstration that higher-order cardioids emerge in the special case of maps with Z c i-i symmetry: in this case 
the symmetry requires that 7 P = and l|5.6p gets substituted by a more sophisticated expression p0.2p . investigated 
in sfTU] (emerging shapes are somewhat less ideal than for d = 2, deviations can reach tens of percents). 
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Figure 30: Mandelbrot Set from Fig[T]with arrows, pointing at particular elementary domains, represented in the table in s|6] Trails 
are also seen in pictures with increased resolutions. Trails are densely populated by clusters. The theory of trails remains an open 
subject. An interesting task is to study how a trail between some two clusters is formed, when we travel across Universal Mandelbrot 
Space as in s l4.ll and these two clusters emerge from a splitting of a single cluster. Another interesting problem is to find the locations 
of the crossing points, where many trails meet together, examples of such crossings are points A and B. 



5.3 Comments 

Note that at the third step we kept terms up to a 2 in expansions of the c-dependent functions, like we did at 
the first step with the x functions. As usual for this type of method to work successfully it is important to 
correlate all approximations: attempt to make one part of calculation more accurate than another decreases the 
total accuracy. Note also that keeping linear terms only, without quadratic corrections, would make eq. (|5.2p 
senseless, and according to the just formulated mnemonical rule one is forced to keep a 2 terms as well. And 
indeed, neglecting them would provide a disaster in description of descendant domains: we already learned in s !2.5 
that their characteristic difference from the root domains is that F p (x, c) — f p (c) + x 2 j(c) + 0(x A ) should vanish 
somewhere at the boundary, while in neglect of the a 2 terms this would be very difficult to achieve while keeping 
fp = fp( c p) 7^ in the center c p of the domain. In fact, the difference between descendant and root domains 
is exactly in <7 2 -terms: their relative magnitude is measured by parameter £ p , and it is negligibly small for root 
domains and close to unity for descendant ones. 



6 On accuracy of the small-size approximation for the family / = x 2 + c 



Numerical characteristics of the lowest (in divisor forest) elementary domains of the Mandelbrot set M.2 are 
represented in the following table (positions of these domains are marked by arrows in Fig |5U|) . 
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V 


1 


2 


3 


3 


4 


4 


4 


4 


c p 





-1 


-1.754877667 


-0.1225611669 
±0.7448617670i 


-1.940799807 


-0.1565201668 
±1.032247109i 


0.2822713908 
±0.5300606176i 


-1.310702641 


distance 
from the root 




(1) 


1 

(2,1) 




(3) 


1 

(3±,1) 




(4i) 




(4a±) 


1 

(4±,1) 


2 

(4,2,1) 


fp 


1 


-1 


-5.649435914 


-1.67528205 
Tl. 1245590i 


-25.53361247 


-9.826826127 
=Fl.391722418i 


-2.273407347 
=F2.878229429i 


1.734079638 


fp 





2 


17.89661552 


-5.9483077 
±6.7473541i 


247.9985718 


-13.82424296 
±86.41684550i 


-30.39333448 
=F9.578562508i 


-15.56341649 




1 


-2 


-9.29887185 


-1.35056409 
T2.249118i 


-39.49472178 


-10.55453437 
=p5.448066568i 


-0.142465954 
=F3.098932717i 


4.88872230 


% 





2 


22.91612617 


-7.458063 
±3.767907i 


-39.49472178 


-42.68642404 
±80.90620560i 


-14.72856928 
=Fl6.79943562i 


-21.82510289 


£ 

from 115.811 


n 
u 

< i 


"I 
1 

= 1 


O 0770^901 1 HQ 

< i 


1 0919^1^ 

±0.04763015i 

« 1 


o oi i z&7f\f\7C\f{'z 

U.UloD i UU I UOo 

< i 


O 0^849^^0749 

±0.08449808566i 

KpI « i 


±0.07045065295i 
w 1 


1 0^^0(^7971 
Si 1 


2r P = (/p7p) _1 
from l|5.7|) 


i 


0.5 


0.0190355 


0.009518 
T0.18867i 


0.00099163 


0.0069178 
T0.0049095i 


0.066394 
T0.057585i 


0.1179602 


c p+ 


0.25 


-0.75 


-1.75 


-0.125 
±0.6495190528i 
( — vx^Fii /r\ 


-1.940550789 


-0.1547246055 
±1.031047228i 


0.25 
±0.5i 


-1.25 


Cp — 


-0.75 


-1.25 


-1.768529153 


-0.1157354238 
±0.8379990280i 


-1.941537753 


-0.1613575037 
±1.036031085i 


0.3161758500 
±0.5574717760i 


-1.368098940 


Cp-j- Cp — 


1 


0.5 


0.018529153 


-0.0092645762 
=F0.1884799750i 


0.000986964 


0.0066328982 
T0.004983857i 


-0.0661758500 
=F0.057471776i 


0.118098940 


(^p+ Cp^K 


1 


2 
0.5 


4 

0.019510668 


2 

-0.0048776662 
=F0.1906854280i 


0.000996072 


0.0071822452 
T0.004799524i 


2 

-0.0645427816 
T0.0601212352i 


2 

0.121405282 



The entries in the last two rows should be both compared with 2r p , which is calculated within SSA in the middle 
part of the table. Since the shapes of root and descendant domains are different, parameters n in the last row are 
also different: n — 4 for cardioid-shape root domains and K = 2 for circle-shape descendant domains. 
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V 


5 


5 


5 


5 


5 


5 


Cp 


-1.625413725 


-1.860782522 


-1.985424253 


0.3592592248 
±0.6425137371i 


-0.04421235770 
±0.9865809763i 


-0.1980420994 
±1.100269537i 


distcLUCG 

from the root 




(50 




(5a) 




(53) 




(54±) 

[...,8,4 ± ,1] 




(5s±) 
[...,12,6,3±,1] 




(5a±) 
[...,12,6,3±,1] 


f P 


-12.346786 


27.952811 


-106.51134 


-15.264582 
=F4.049590i 


-1.4606030 
Tl5.65788i 


-34.451623 
±8.245915i 


fp 


32.445339 


-211.46971 


3811.7679 


-270.4039 
±144.8293i 


-190.8108 
=Fl51.4911i 


421.6135 
±736.0759i 




-19.95443 


45.84473 


-161.34688 


-9.4954135 
=F9. 112127i 


6.050675 
Tl7.26642i 


-41.054355 
=F3.02090i 


> 


17.8489 


-272.6037 


5625.2834 


-237.5248 
T0.045178i 


-108.0265 
T230.0270i 


207.2558 
±911.8688i 


from 15.81 


u . uuoy doo^o 

< i 


n 007^0.1 Szin 
u.uu i oy io^tu 

« i 


u.uuouoyy < 

< i 


u.uzozoyoz 
±0. 1305419* 
I€pI « 1 


U.UOoOO I oo 

=F0.06450516i 
\(p\ « 1 


n nn^i i 79s 
±0.02358242i 
\b\ « 1 


2r p = (/p7p) _1 
from 115.711 


0.004058884 


0.0007803423 


0.0000581894 


0.00250125 
=F0.00411026i 


0.0033726 
±0.00083981i 


0.00067682 
±0.00011025i 


c p+ 


-1.624396989 


-1.860586973 


-1.985409691 


0.3599331332 
±0.6415066668i 


-0.04506136598 
±0.9868115622i 


-0.1978729467 
±1.100298438i 




0.004066944 


0.000782196 


0.000058248 


0.0026956336 
T0.0040282812i 


0.00339603312 
±0.0009223436i 


0.0006766108 
±0.000115604i 
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p 


5 


5 


5 


6 


6 


6 


6 


Cp 


-1.256367930 
±0.3803209635* 


0.3795135880 
±0.3349323056* 


-0.5043401754 
±0.5627657615* 


-1.476014643 


-1.907280091 


-1.966773216 


-1.996376138 


distance 
from the root 




(5r±) 

[. ..,12,6±,2,1] 


1 

(5i±,l) 


1 

(5 2 ±,1) 




(6i) 




(62) 




(63) 




(64) 


fp 


-9.720195 
Tll-81330i 


-2.700383 
=F5. 404227*' 


-3.449207 
Tl. 112266* 


9.557119 


-73.91417 


135.0997 


-431.94389 


fp 


-59.44225 
±171.0187* 


-48.15074 
T78. 38748* 


9.823111 
±30.87527* 


-106.133658 


491.1284 


-2562.839 


60341.247 


7 P 


-5.71067 
^21.95456* 


0.809902 
=F3. 399590* 


-2.871463 
T2.098768* 


20.23824 


-115.0273 


207.856 


-649.8852 


jp 


-165.8458 
±183.4574* 


2.830436 
=F47.49326i 


0.04567 
±29.22689* 


-150.6058 


576.3060 


-3609.262 


90085.134 


from H5.8H 


0.0176798 
=F0. 0451181* 
I?pI « 1 


1.00083 
±0.0866142* 
»s 1 


0.984283 
T0.000539978* 
« 1 


0.067182 
< 1 


0.00098004 
< 1 


0.00095613 
< 1 


0.0007426 
< 1 


2r P = (/p7p) _1 
from 115.711 


0.001692541 
=p0.00233202i 


0.04612246 
T0.0107757i 


0.04556086 
=F0.0627926i 


0.005170116 


0.000117617 


0.000035611 


0.0000035623 


c p+ 


-1.256801994 
±0.3797412022* 


0.3567627458 
±0.3285819450* 


-0.4817627458 
±0.5316567552* 


-1.47469537780 


-1.90725067795 


-1.96676431090 


-1.99637524690 


K 

(Cp+ - Cp)K 


4 

0.0017363 
^0.00231905i 


2 

-0.04550168 
=F0.01270072i 


2 

0.04515486 
=F0.06221801i 


4 

0.00527706 


4 

0.000117652 


4 

0.000035620 


4 

0.000003564 



4G 



p 


6 


6 


6 


6 


6 


6 


Cp 


0.4433256334 
±0.3729624167i 


0.3965345700 
±0.6041818105i 


0.3598927390 
±0.6847620202i 


-0.01557038602 
±1.020497366i 


-0.1635982616 
±1.097780643i 


-0.2175267470 
±1.114454266i 


distelllCG 

from the root 




(65±) 

[...,10,5i±,l] 




(6b±) 
[...,8,4±,1] 




(6r±) 
[...,8,4±,1] 




(6s±) 
[...,6,3±,1] 




(6 9 ±) 
[...,6,3±,1] 




(6io±) 
[...,6,3±,1] 


f P 


-22.131316 
T8.549589i 


-7.730347 
T29.76415i 


-45.44965 
±11.11460i 


-47.08247 
=p42.3870i 


-13.70794 
T59.0423i 


-94.05752 
±65.91697i 


f P 


-868.64319 
T90.27354i 


-427.073 
Tl083.050i 


-431.8717 
±2162.074i 


-3003.383 
±437.063i 


-2600.199 
=F913.171i 


7694.999 
±2405.039i 




-8.156137 
=pll.41327i 


6.87116 
=F23.4912i 


-37.7109 
T9.24215i 


-28.2862 
=F64.5198i 


7.68865 
=F69.5663i 


-125.4789 
±40.9618i 


ip 


-389.998 
=F307.274i 


151.331 
=F866.959i 


-1097.368 
±1391.979i 


-3129.19 
=F931.571i 


-2370.323 
Tl950.970i 


7560.924 
±5049.389i 


t 

?p 

from II5.8II 


n ooao^ 
±0.16159i 
I€pI « i 


0791 ^ 

T0.01058i 

l«pl « 1 


o 01 f;fis 
±0.034629i 

I£pI « i 


±0.005837i 
I«pI « 1 


n ni 0090 
=F0.01208i 
\b\ « 1 


O OOI 788 

±0.006621i 

\Zp\ « i 


2r P = (/p7p) _1 
from 115. 71 


0.0007487 
=F0.0029099i 


0.0013280 
±0.00004046i 


0.00055046 
T0.000000276i 


-0.00007044 
T0.0002127i 


-0.00023408 
=F0.00002776i 


0.000039602 
±0.00005275i 




0.44355069141 
±0.37224682H 


0.39619446067 
±0.60419336859i 


0.3600296164 
±0.684763498* 


-0.01558797607 
±1.0204438975i 


-0.163657003678 
±1.0977739136i 


-0.217516881 
±1.114467467i 




0.000900232 
=p0.00286238i 


2 

0.001360437 
±0.00004623i 


2 

0.00054751 
±0.00000591i 


0.000070360 
=F0.00021387i 


-0.000234968 
=F0.00002692i 


0.000039463 
±0.00005280i 
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V 


6 


6 


6 


6 


6 


6 


Cp 


-0.5968916446 
±0.6629807446i 


-1.284084926 
±0.4272688960i 


0.3890068406 
±0.2158506509i 


-1.772892903 


-0.1134186559 
±0.8605694725i 


-1.138000667 
±0.2403324013i 


disteincG 
from the root 




(6u±) 
[...,10,5 2± ,1] 




(6l2±) 
[..., 12,6±,2, 1] 


1 

(6±,1) 


1 

(6,3) 


2 

(6,3±,1) 


2 

(6±,2,1) 


f P 


-21.17430 
Tl.239334i 


-49.4808 
=F43.51796i 


-2.84507 
=F8.81348i 


6.02446709 


2.710416 
±2.01496i 


3.0373845 
±0.674370i 


fp 


258.9018 
±390.3375i 


-401.263 
±2355. OOli 


-13.2205 
=F228.2727i 


-699.2133940 


53.81286 
T56.37906i 


-9.154561 
=F50.431515i 




-19.86933 
=F7.65439i 


-40.9497 
=F85.4133i 


1.51341 
=F3.51012i 


19.09634167 


3.35657 
±5.72766i 


3.882374 
±4.659455i 


> 


122.1821 
±435.7183i 


-1796.0107 
±2809.8473i 


47.3746 
=F74.1184i 


-1106.4705637 


71.1246 
=F30.8637i 


28.408033 
T64.35131i 


ip 

from 115.811 


U.UDZDD I 

±0. 034437i 
\b\ « 1 


U.UUO I ( O 

T0.006297i 

I€pI « i 


±0.09882i 
1 


l.UUoUD 

w 1 


±0.013004« 

ft! 1 


±0.0438« 

ft 1 


2r P = (/p7p) _1 
from 115, 71 


0.0020161 
=F0.0009153i 


0.000043399 
=F0.00015422i 


0.0281207 
±0.0026745i 


0.008692230 


-0.0048602 
=F0.044335i 


0.0242924 
=F0.047098i 


Cp+ 


-0.5963742242 
±0.6627532528i 


-1.284095877 
±0.4272302898t 


0.375 
±0.2165063509* 


-1.768529153 


-0.1157354238 
±0.837999027i 


-1.125 
±0.2165063509i 


( C p+ Cp)K 


0.00206968 
=F0.000909967i 


0.00004380 
=p0.00015442i 


2 

-0.028013681 
±0.001311400i 


2 

0.008727500 


2 

0.0046335357 
=F0.045140890254i 


2 

0.0260013340 
T0.0476521007i 



A few comments to the table are now in order: 

We present the values of parameters with high accuracy, which strongly exceeds our needs in the present paper, 
but this data can be used in the future investigations. One should not be surprised by the high accuracy of 
experimental data: since it is computer experiment over Platonian entity, accuracy is unlimited. Moreover the 
numbers in the last column can be also reproduced as resultants zeroes [2]: though it is a difficult calculation 
(beyond capabilities of MAPLE on an ordinary laptop already for p > 6), its accuracy is in principle unlimited. 

c + is closer to the root than c_ : c + is the point of merging with the parent domain or the cusp position if the 
domain is itself a root, while c_ is the "opposite" point, i.e. the merging point with the next descendant of the 
order 2. 

Starting from p = 5 there are many root domains of the type (5), (6) etc, and - starting from p = 3 - many 
descendants (3, 1), (4, 1),. . . , (6, 2, 1) etc: a-parameters begin to emerge. The two domains (3±, 1) differ by complex 
conjugation only, but in other cases systematization in a-sector is less straightforward, their sizes and orientations 
depend essentially on a. Still, because of the symmetry of the Mandelbrot Set under complex conjugation, the 
domains with centers at non-real c come in pairs. Such complex conjugate domains are always labeled by indices 
±. 

Positions of the domains in divisor forest are shown in the second line of the third row. For root domains the 
original direction of the trail, connecting it to the central cluster, is also shown in square brackets in the third row. 
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Of course, all root domains with centers at real values of c belong to the trail, originating at [. . . , 2™, . . . , 16, 8, 4, 2, 1], 
and it is not mentioned in the table. 



From this table we observe: 

- the good accuracy of the relation (see the last three columns) 

2r p ~ ^p-f Cp— 

between the theoretically-predicted (in the small-size approximation) complex-valued size r p of an elementary 
domain and the difference between experimentally found extreme points c p+ and c p „ ; 

- the correlation between the value of £ p and the distance of elementary domain from the root of the corre- 
sponding cluster (the corresponding columns are boldfaced): £ p is tiny for the roots (distance = 0) and close to 
unity for all descendants (distance > 1). 



7 Why £ p « 1 for descendants: la raison d'etre for circles 
7.1 Approach to description of descendants 

We can study descendants of a given elementary domain within the same small-size approximation (SSA), simply 
iterating approximate expression 

f op ^)^f P (c)+lp(c)x 2 

to 

f° {2p) {x) « f 2p {c) + l2p (c)x 2 « f p (c) + 7p (c)(/p(c) + lp {c)x 2 ) 2 » f p (c)(l + f P (c) lp (c)) + 2f p (c) 1 2 (c)x 2 (7.1) 
and so on. Thus in this framework 

72 P (c) « 2/ p (c) 7p 2 (c) (7.2) 

C2 P is a non-trivial root of this new /2 P (c), 

fp(c 2 php(c 2p ) = -1 (7.3) 

This procedure - if at all justifiable - can be valid only for C2 P , associated with a descendant domain of c p (but 
not a root domain of some new cluster), since it relies on SSA and assumes that C2 P is very close to c p . The shift 
a 2p = C2p — c p can actually be found in SSA by solving (17. 3p iteratively: 

"Wl-^^-l (7.4) 



2t* p y 27* p j 

Now we are going to demonstrate that evaluated for such c 2p within SSA, is indeed equal to unity (this is 
no more than a consistency check, because validity of the SSA itself will not be theoretically justified). Afterwards 
this calculation is extended to descendant c mp for all to. Further, eq. (|7.4p and its generalizations for a mp are used 
in s[5]to evaluate SSA approximations of various Feigenbaum indices. Finally, in s llOl we briefly consider the case 
of specific ^d-i-symmetric f(x; c) — x d + c families. 



7.2 Evaluation of £ 2p for a descendant 

This is a rather straightforward calculation. From (|7.1[) we obtain - in the small-size approximation, after substi- 
tution of c = C2p and (|7.3[) . and after expanding functions of C2 P in powers of <72 P = C2 P — c p from (|7.4[) - the set of 
recurrent expressions: 



/ 2 p(c) = / p (c)(l + 2/ p (c) 7p (c)) + f 2 (c)j p (c) 



-fp(C2p) ~ fp{c 2 p) — {c 2 p) 

7p 




"fv 




) 


jp 


7p_ 





(7.5) 
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72 P (c) = 2/ p (c) 7p (c) W 72 P = 72 P (c 2p ) = -27 p (c 2 p) = -2 7p + cr 2 p^J (7.6) 

/ 2p (c) = / P (c) (l + 2/ p (c) 7p (c)) + 2/ p 2 (c) 7p (c) + 4/ p (c)/ p (c) 7p ( C ) + %(c)%(c) =j. / 2p = 2/ 2 7p (7.7) 

7 2p (c)=2/ p (c)7 p 2 (c)+4/ p (c)7 p (c)7p( C ) => 72 P = 2/ p (c) 7p 2 (c) (7.8) 
Substituting these expressions into (|5.7p and (|5.8[) . we obtain: 



and 



as required. 



,2p 



11 r 

r 2p = — i ~ p 77 7T = — 7 ~ T (T.9) 

2/ 2p7 2p 4 / p7p {1 + ^(1+2^)} 2(1-^) 



hp , 72p^ _ 1 / 2/ p 2 7p | 2/ p (c) 7 2 (c)^ _ 2 / p7p _ ^ ^ io ^ 



/2p 7 2p \2/2p 72p / / 2 p 7 2p \ — 2/p 2 7p / /2p72p 



7.3 The rules of SSA 

Note, that within SSA we consider r p f p / f p and r P 7 P /7 P as small parameters and ignore their quadratic powers as 
well as higher derivatives. This is needed for self-consistency of the SSA, even despite individual corrections need 
not be small (especially for domains which are not the first descendants, i.e. when £ p w 1 is not small) - however, 
if included, they should come together with other corrections to the SSA, which were also ignored. Actually, as we 
saw in sJS] the summary effect of all corrections is small, but the theoretical reason for this conspiracy in the case 
of higher descendants remains to be identified. 

It deserves formulating the rules of SSA explicitly: 

• Expand in powers of x and leave the first two non-trivial terms (constant and x 2 in the case of x 2 + c family) 
- for generic value of c. 

• Expand in powers of a — c — c cr n and leave only the first corrections ~ /// and ~ 7/7. 

• If two different but two close c cr ;t appear in the problem (say, centers of two adjacent elementary domains), 
expand in powers of their difference, leaving only the first two powers of the difference. 

• Combining all these expansions, keep only the first two corrections in expressions for the final quantities, in 
practice this means keeping all powers of / and 7 and ignore everything beyond the first powers of /// and 7/7. 

7.4 Position and radius of arbitrary descendant domain 

Generic descendant domain has a parent of order p and has itself a multiple order mp. It is attached to the parent 
at a zero of the resultant R x (^F mp (x)/ F p (x), F p (x)j . Parent can be itself a descendant and a chain of ancestors 
lead to a root domain of the cluster, however, only the first term in this chain - the mother domain, of which the 
domain of interest is an immediate descendant, - is relevant in the SSA-based calculations. 
It is easy to check that generalization of the SSA relation (17. 2ft to arbitrary m is 



fmp(c) ~ 7-^/m(/p7p(c)) (7.11) 

7mp(c) ~ 7p(c) 7m(/ P 7p(c)) ( 7 - 12 ) 



Now descendant root c mPjP of f mp is defined by the choice of immediate descendant c m i for f m : from 

/™(c m ,i)=0 (7.13) 

we have 
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For comparison with s !7.2l one should keep in mind that for / = x 2 + c there is a single order-two critical point 
C2,l = c 2 = -1. 

Note, that not all the zeroes c m of (|7.13p describe immediate descendants (m, 1) of the central domain (1): some 
provide the new root domains (m) or higher descendants (m, mi, ■ . .) with nontrivial divisors mi ^ 1 of m. These 
extra zeroes (especially associated with domains from the different clusters) should not be used in the following 
calculations, because they correspond to remote domains and SSA has no reason to work for them. 

Repeating for generic m the calculations, performed s !7.2l for particular case of m = 2, we obtain: 



fm P (c) ~ ^^y/m(/ P 7p(c)) (/ P 7p(c) + / P 7p(c)) - ^f|^y/™(/p7 P (c)) (7-15) 
In combination with (|7.12[) this implies that 

^ — fmp r Ymp(Cmp,'p) ~ ,/*m7m(Cm,l) ( fp^fp ~t~ /p7p ) (Cmp,p) ~ /m7m(^m,l) ( /p7p ~t~ ^rap.p ( /p7p ^fp^p ) ) (^p) 

In the first transformation we omitted one term with f m {c m ,i) — 0, and in the second transformation we defined 
functions at c mP)P through their values at c p , keeping only the first non-trivial term of Taylor expansion in powers 
of a m p,p = Cmp.p — c p . This shift is defined in a similar manner from (|7. 14|) : 



Vmp^fplpiCp) + 2 fJ mp,p \ fpTp fplp^j ( C p) ~ c m,l 



(we remind that f P {c p ) = 0). Substituting for remaining parameters f m j m — (2r m ) , f p j p — (2r p ) 1 and 
/p7p + 2/ P 7p = — ^p(2rp) _1 , we obtain for the counterparts of (|7.4[) and (|7.9p : 

and 

r mp » 2r m r p (\ - £ p ^^) (7.17) 

7.5 Evaluation of generic £ mPiP 
For evaluation of £ m p,p we need also 

fmp(c) « -I-^/ m ^/ p7p (c)^ ((/p7p)) + ^^y/m(/p7p(c))'('/p7p) ~ 

-2^4/ m (/p 7 p(c))(/p7p) + ^\fJfplp{c)) + 0(7) (7.18) 
7p \ c ) 1 P \ c ) 

and 

7m P (c) ~ 7 P (c)7m (/ P 7p(c)) (/p7p(c) + /p7p(c)) + 7p(c)7m (/p7p(c)) (7.19) 
At c = c m p jP the terms with / m (c m ,i) = do not contribute, and we obtain 

/l ■■ • \ /l •• • . \ / • \ 2 

) ~ ( 7m + /m7m J (C m ,l)(^/p7p + /p7pj (c mP: p) + 

+ fm7m(c m ,l) Q/p7p + /p7p^ {c m p,p) (7.20) 

In order to get £ m p,p we divide by the square of 

fmp^fmp ~ fm'Yrn(Cm 1 l) ^/p7p ~t~ /p7p^ (^mp,p) (^"21) 
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and change sign, so that (|7.10p generalizes to and 




In 



fmplmp \ 2fmp Imp J fmlm \ 2fm 7m J 

fp % \ 1 t 2r m £ p mjtl 



fmHm \2/p 7p / (/pTp fp'jp) ipmp) 



2 



Cm (7.22) 



Keeping the second term at the r.h.s. is beyond the accuracy of the SSA and it should be neglected (we ignored 
it in (|7.10p . but kept in (|7.22p to preserve formal consistency with the case m = 1, when r\ = |, £1 = and, of 
course, a PtP = 0). 

From (|7.22[) it is clear that if £ mj i « 1 for direct descendant (m, 1) of order m of the central root domain 
(ci = 0), then in the SSA £ mp , p ~ 1 for all other descendants, at all levels in all clusters. In s !7.2l we exploited the 
fact that for m = 2 the r.h.s. is extremely simple: for /2(c) = c(c+ 1), 72(c) = 2c and C2 = —1 it is obviously unity. 
In s[6]we saw that £, m ,i is indeed close to unity for m < 6, and it is natural to believe that this remains true for all 
m, however no theoretical explanation of this fact is yet available. Still, if accepted, it implies that £mp,p ~ 1 for 
all p > 1. 



8 Feigenbaum indices 

8.1 The case of period-doubling, m = 2 

It is now time to solve quadratic equation (|7.9|) : the distance between the centers of a parent domain (p, . . .) and 
its immediate descendant (2p,p, . . .) is 

{— 2r p if £ p « i.e. for the first descendant 
(8.1) 
(1 — \ // 5)r p if ps 1 i.e. for a higher descendant 

Substituting this into (|7.9p . we obtain: the radius of descendant domain (2p,p, . . .) is 

!2?"p if £ p « i.e. for the first descendant 
if £p ~ 1 i.e. for a higher descendant 
Thus we get for the Feigenbaum doubling parameter 82 = Virci p ^r X1 (r p /r2 P ) 

8 2 « 2V5 = 4.4721 .. . (8.3) 

(exact value is known to be 82 = 4.6692 . . .). Note that 82 approximately acquires this value already for the second 
descendant of the root, far before the p — > 00 limit. 
Consistency requires that 

C2p + r 2p = c p - r p (8.4) 

i.e. 



o-2p = r 2 p + r p w r v [ 1 + ^= ) (8.5) 



This is indeed almost true: 1 + ^ w V5 - 1 (1.2236 . . . « 1 



2361. 



The west-limit point Coo of the central cluster (superscript is (2) because the point is obtained by a sequence 
of doublings of the order of the orbits) can be represented as 

c£»-\-2r2-2r,-... = -\-2Yr2 km = -\--^ = -l.mm... (8.6) 

k=l 2\/5 

(we remind that for the first descendant r 2 = \r\ = 4) or, alternatively, as 



eg) = c 2 + £ a 2k+1 « c 2 - - -1 - = 3980 ■ ■ • 

fc=i 1 - 275 2 2V5 - 1 

The difference between these two values characterize accuracy of the SSA, and within such error they coincide with 
exact value c^, = —1.4012 . . .. 
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8.2 The general case (arbitrary m) 

Solving ()7.16|) . we obtain: 

2c m r p if £ p w i.e. for the first descendant 



' rap.p 



' rap 



(1 — yl — 4c m )r p if ^ w 1 i.e. for a higher descendant 

^TmXp if £p ~ i.e. for the first descendant 
Ji-Ic ^ ~ 1 i- e - f° r a higher descendant 



(8.9) 



Thus we obtain for the Feigenbaum parameter S m = limp—^ (r p /r mp ) 

VI - 4c m 
m ~ 2r 

and for the complex-valued ratio £ m = <r m p p/rp we get 



(8.10) 



e m = 1 - VI - 4c m (8.11) 
When (p, . . .) is itself a descendant domain and has circle rather than cardioid shape, the consistency condition 

\<J mp \ = Vp\ + \r mp \, (8-12) 
expressing the distance between centers of two touching circles through their radiuses, implies that 

\e m \ = l + \S m r 1 (8.13) 



or 



2|r„ 



|1-VT^4^| = 1 + J-^—. (8.14) 

A more detailed consistency condition includes not only distances, like (|8.12|) . but also exact position (the phase 
<p m ) of the touching point between the circles (mp,p, . . .) and (p, . . .): 



mp ~r~ r nl p c- / p ~r ^p 



i(f> 



e l(pm r p + c p (8.15) 



This means that 



e m = e**» - ^- (8.16) 
The end-point of an infinite sequence of descendant domains (p), [rnp, p), (m 2 p,mp,p), ... is given by 



Mm) _ a k - c + f W * -r I £mr? "P - r + 2r lr I r? " £m 1 ^8 17^1 

°oo °p i / u m^p b mp ~ c m / 1 m K p u mp i 1 c— 1 P P I m -p r— 1 J \ / 

fc=i fe=i l - d TO v i - o w 



In particular, for the central cluster with (p) = (1) 

C S m) = + T^T (8-18) 

Within SSA the only input in all these formulas for a given m consists of two complex numbers: c„, = c m ,i 
and r m = r m .i, characterizing the properties of the next-to-root domain (m, 1) in the central cluster. These c m 
and r m are entries of the table in s[51 Taking c m and r m from that table, we now make a new one, comparing 
predictions of eas. ([8T8 |) - ((87T8]) with experimental data. Numbers in square brackets in the last column are positions 
of the limiting points, measured with the help of Fractal Explorer. 
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domain 
(m, 1) 


Cm 


1r 
Zl ' m 


e i<t>m 




(fmp / Tp) 


<5m 
VP 1 ^mp ) 






•-'GO 

from J8.18J 


\J 1 - 4c m 




l|8.16|l 


(2,1) 


-i 


0.5 


-1 




1 - v 7 ! 




VE- 1 

1.2361 w 1.2236 


l-VE 


I 15— y/E 

= -1.3980 
[-1.401] 


(3,1) 


-0.123 
+0.745? 


-0.01 
-0.19? = 
0.19 • e 1 ' 03 ^ 


-0.5 
+0.87? 


1.55 
-0.96? 


-0.55 
+0.96? 


4.61 
+8.42? = 

9.59 ■ e 102 ^ 


1.107 
S3 1.104 


-0.553 + 0.959? 
w -0.550 + 0.957? 


-0.020 
+0.785? 
[-0.0234 + 0.7836i] 


(4,1) 


+0.282 
+0.530? 


-0.066 
-0.06? = 
0.089 • e - 102 ^- 


i 

04 - — 


0.999 
-1.061? 


0.001 
+ 1.061i 


-0.2848 
+16.34? = 
16.34- e 1 ' 01 ^ 


1.061 
sa 1.061 


0.001 + 1.06? 
pa 0.001 + 1.06? 


0.3115 
+0.4932i 

[0.3098 + 0.4947?] 


(5i,l) 


0.380 
+0.335? 


-0.046 
-0.011? = 
0.047 ■ e -°- 937ri 


0.309 
+0.951? 

<PBl = ~ 


0.677 
-0.989? 


0.323 
+0.989? 


-9.06 
+23. 7? = 
25.35 ■ e 3 08 "^ 


1.041 
w 1.039 


0.323 + 0.989? 
ss 0.323 + 0.988? 


0.377 
+0.311? 
[0.3770 + 0.3117i] 


(52,1) 


-0.504 
+0.563? 


0.046 
-0.063? = 

0.078 ■ e -° m7 ^ 


-0.809 
+0.588? 

<t>5 2 = — 


1.841 
-0.612? 


-0.841 
+0.612? 


20.25 
+14.44? = 

24.87 -e ^ 


1.040 
sa 1.040 


-0.841 + 0.612? 
-0.842 + 0.611i 


-0.503 
+0.605? 

[-0.5031 + 0.6048?] 


(6,1) 


0.389 
+0.217* 


-0.028 
+0.003? = 
0.028 ■ e-°- 977ri 


0.5 
+0.866? 

06 — 


0.487 
-0.891? 


0.513 
+0.891? 


-20.57 
+29.61? = 
36.05 • e 2 ' 08 TT 


1.028 
ss 1.028 


0.513 + 0.891i 
S3 0.516 + 0.889? 


0.380 
+0.206? 
[0.3810 + 0.2047?] 



Of course, one can consider limiting points of other sequences, not obligatory of the type [. . . ,m,m, m]. One of 
the open questions is if there is any difference between periodic (after some step) and aperiodic, i.e. "rational" and 
"irrational" sequences. Another important class consists of sequences [. . . , 2, 2, . . . , 2, m r , . . . , mi], ending by 2's 
only - they describe normals to the cluster's boundary and serve as origins of trails, connecting the cluster with 
its neighbors. 



9 Cardioids and resultant zeroes 

As explained in [2J, a boundary of domain (p, . . .) is densely populated with a countable set of its merging points 
with descendant domains (mp,p, . . .), located at zeroes of the resultants R{G mp ,G p ) with all integer m. Since 
within SSA the boundaries are well approximated by cardioids and circles, and merging points are characterized 
by the angles <f> m = — , one can expect that simple approximations exist for locations of the resultant zeroes in 
terms of c p , r p and e ik ^ m with k = 1, ... ,m — 1. This is indeed the case, at least for the Mandelbrot Set, i.e. the 
family {f(x)=x 2 + c}. 

For example, the zeroes of R(G m , G\) — they can be found among the values of c p + in tables in s(6] — are given 

by 

The first values of this quantity are: 





(2) 


(3) 


(4) 


(5i) 


(52) 


(6) 


(7i) 


(72) 


(73 ) 


c (1) 


-0.75 


-0.1249999999 
+0.6495190530? 


0.25 
+0.5? 


0.3567627456 
0.3285819454? 


-0.4817627458 
0.5316567550i 


0.375 
+0.2165063510? 


0.3673751344 
+0.1471837632? 


0.1139817500 
+0.5959348910? 


-0.6063568845 
+0.4123997402? 
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When k are not shown, it is equal to unity, k = 1. Similarly, the zeroes of R(G2m, G2), belonging to the boundary 
of descendant domain (2, 1), are given by: 

« c 2 + r 2 e 2mk ' m = -1 + (9.20) 





(2) 


(3) 


(4) 


(5i) 


(5a) 


(6) 


(7i) 


(7a) 


(7 3 ) 




-1.25 


-1.125 
+0.2165063510i 


-1 
+0.25« 


-0.9227457516 
+0.2377641291i 


-1.202254249 
+0.1469463130i 


-0.875 
+0.2165063510i 


-0.8441275496 
+0.1954578706i 


-1.055630233 
+0.2437319780i 


-1.225242217 
+0.1084709348i 



and zeroes of -R(G 3m , G 3 ), belonging to the boundary of descendant domain (3±, 1)) - by 

Cmf 1] ~ c 3 + r 3 e ±2Trlk/m = -0.1226 ± 0.7449? - (0.0047 ± 0.0943z)e ±2T " fe/m (9.21) 





(2) 


(3) 


(4) 


(5,1) 


(5,2) 


(6) 


(7,1) 


(7,2) 


(7,3) 


(3±,l) 


-0.118 
±0.839i 


-0.039 
±0.788i 


-0.028 
±0.740i 


-0.034 

±o.7iii 


-0.063 
±0.818i 


-0.043 
±0.694i 


-0.052 
±0.682i 


-0.030 
±0.761i 


-0.077 
±0.828i 



while those belonging to the boundary of the root domain (3) are 

/ 1-Kikjrn \ / 2-Kik/m \ 

~ c 3 + r 3 e 2mk /" 1 I 1 - ^— J = -1.7549 + 0.0095e 2mfc/m I 1 - =— J (9.22) 





(2) 


(3) 


(4) 


(5i) 


(52) 


(6) 


(7i) 


(72) 


(7s ) 


c (3) 


-1.769 


-1.757 
+0.012i 


-1.750 
+0.009i 


-1.748 
+0.006i 


-1.764 
+0.010i 


-1.748 
+0.004i 


-1.748 
+0.003i 


-1.753 
+0.011i 


-1.767 
+0.008i 



Since domains (3, 1) and (3) are circle and cardioid only approximately, accuracy in the last two tables is relatively 
low and we do not keep as many digits as in the first two tables. Still, the numbers in the tables reproduce actual 
positions of resultant zeroes at percent-level accuracy, standard for the SSA in the case of the Mandelbrot Set. 
Thus, not only the shapes of elementary domains are nicely represented by cardioids and circles, but all the merging 
points of stable orbits at the boundaries (zeroes of the corresponding resultants, [2]) can be easily found by the 
SSA methods. 



10 The case of Z^-symmetric maps f(x; c) = x d + c 
10.1 SSA in the case of Zj i-symmetry 

In this case all the iterated maps are expanded in powers of x d and in SSA we truncate them as follows: f op (x; c) — 
fp{ c ) + x dr y p (c) + 0(x 2d ). Then the boundary of elementary domain, surrounding a root c p of f p (c), is defined by 

(/ p (c)=x(l-x« 7p (c)) (1Q1) 
\ d-x d - 1 lp {c)=e^ d - 1 ^ 

or, as generalization of (|5.3p . 

p »(d-i 



f P {c)(d lp {c)) d 1 = e<+ (l - e — r - J (10.2) 
Now we need to expand the l.h.s. in powers of a = c — c p and leave the first d terms of the expansion. 
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10.2 Example of the p = 2 domains for d = 3 and d — 4 



We consider here the first descendants of the central elementary domain c = in the case of d = 3 and d = 4: 
relations like (17.22|) should be used to extend the result to all other descendants. Also we restrict example to p = 2 
only. 

From f° 2 (x, c) — (x 3 + c) 3 + c we read: 

/a(c) =c(c 2 + l), 

72(c) = 3c 2 (10.3) 

and the critical values C2 = ±i. Eq. (|f 0.2[) now states: 

3c 2 (c 2 + 1) = u(l - ^u 2 ) (10.4) 

with u — e 1 ^, we need to substitute c = ±i + a and check that (|10.4[) is approximately - modulo terms ~ 0(u 4 ) - 
solved by 



a = r-2U 1 — 



a 

2 



- bu l (10.5) 



with negligibly small a and 6. Substitution of this ansatz into (|10.4[) gives a = -C 1 and 6=^<1. Also from 
the same calculation r 2 — s, and this is in good accordance with reality: the first descendant domain (2, 1) + in 
Fig[2]is bounded by the points c + = l,09z and c_ = 0.77z, so that c + — c_ = 0.32z ps 2r 2 = 0.33i. 
Similarly, for d = 4 we have 

2 4 / 3 c 2 (c 3 + l) =« ^1- ^ (10.6) 
and for c = CJ3 + a, CJ3 = e i7r ^ 3 we obtain 

with a = 2~ 4 / 3 = 0.40, b = 11 • 4" 1 / 3 /9 = 0.77 and c = 4/81 < 1. The biggest "diameter" of this elementary 
domain is 2 „ 2 \y 3 (l + |) s» 0.33, in good agreement with c_ = —1.10, c + = —0.78 for the (2,1) domain in 
FigEl The distance between the two cusps of this deformed cardioid is approximately 0.4 of the biggest diameter, 
what is also in agreement with Figj3] (ordinate of the cusp, which is shown by arrow in the picture, is 0.063, and 
2 • 0.063/0.33 f» 0.4). Since b < 1 in (|10.7[) . the cusps have finite angles, what is not confirmed by FigJT] the true 
value of b is close to unity - the difference 1 — b w 0.23 is inaccuracy of SSA in this example. 



11 Conclusion 

In this paper we calculated the shapes of elementary domains of the Mandelbrot set [5], following the general 
algebro-geometric approach of [2]. We explained the qualitative features of these shapes, found the origin and 
number of cusps, explicitly showed how they change when one Mandelbrot set is deformed into another inside the 
unifying Universal Mandelbrot set. We showed that the nearly ideal cardioid and circle shapes of these domains 
in M2 (Figfljl are nicely described in the small-size approximation, based on truncating the relevant polynomials 
to the first orders in deviations x — x cr and c — c cr from their critical values. It is not a big surprise, but some 
conspiracy is needed - and was indeed found in the behavior of parameter £ p , which is not always small, as one 
could naively expect - to explain the coexistence of different structures: cardioids of different orders. 

We did not give a theoretical justification of the small-size approximation - next-order corrections were not esti- 
mated - instead its percents-order accuracy was demonstrated by comparison of its predictions with the properties 
of the actual Mandelbrot set (measured with the help of the Fractal Explorer [4]). Accuracy is actually much higher 
than one could expect from the over-simplified calculations in [5] , for example the small-size-approximation of the 

ordinary Feigenbaum constant S^ 3 ^ = = 4.4714 ... is much closer to experimental value 62 = 4.6692 . . . than 
82 — 4 + \/3 = 5.7321 ... of rcf.[8j. The systematic approach allows to find all Feigenbaum indices in the same 
way, moreover other characteristics, including continuous, like shapes of elementary domains, not only their sizes, 
are straightforwardly calculated. 
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We demonstrated that characteristics of elementary domains in A^2 are nicely encoded by two parameters like 
r p and £ p , which by recursive formulas like 

r mp w r p ■ (5 _1 (c mj i) (H-l) 

and 

, | for the root domain (112) 
p ~ 1 1 for other domains 

are expressed through the size r root of the root domain in the given cluster and through the critical values c m< i - 
positions of centers of immediate descendants of the central root domain. However, these remaining parameters 
need to be evaluated from sophisticated algebraic equations. As explained in [2], the equations emerge from 
universal structures in particular section of the Universal Mandelbrot set (UMS). Naturally, some characteristics 
of such arbitrary section look arbitrary - at least from its internal perspective. Hopefully, a better understanding 
of r root and c m distributions can be found at the level of UMS, but this remains beyond the scope of the present 
paper. 



It remains to emphasize that investigation of Mandelbrot sets is not just an interesting problem by itself, it 
is crucial for understanding of the future physics, which is going to deal with essentially multi-phase systems, far 
from equilibrium and from the trivial end-points of renormalization groups. One of the main lessons of Mandelbrot 
theory [5] is that phase transitions are not just rare isolated events, concentrated on smooth hypersurfaces in the 
space of coupling constants. Examples of such phase transitions are given by particular merging points between 
two elementary domains (say, between (2, 1) and (1)) - these isolated points in particular A4d in Fig|7] form a 
nice complex-codimension-one hypersurface in UMS (partly represented in Figl27|. However, the true picture - 
FigsflE] - is very different: the entire variety of various phase transitions (mergings of all elementary domains 
of all orders) is not just a collection of particular transition lines. Instead they form a profound new structure, 
moreover they tend to condense and fill entire boundaries of elementary domains, i.e. dimension of the phase 
transitions variety increases as compared to the naive one (and actually its real, not complex codimension in the 
space of complex couplings, is one!). Within particular slices like particular Mandelbrot sets, different phases now 
get fully disconnected, and analytical continuation between them, if at all possible, essentially depends on the 
properties of the new fundamental entity: the UMS, which scientists even did not begin to study! It is the UMS 
that is behind the sophisticated phase structure fTU] of stringy r-functions - effective actions of various multi-phase 
systems, classical or quantum. It is the UMS that one encounters in various problems, from baby-universe creation 
in modern cosmological models to optimization of cooling processes in various solid-state technologies. Still, despite 
its central role in the mysteries of uncertainty, there is no mystery in the UMS itself: it is one of the most important 
and structurized mathematical objects - the universal discrminantal variety, a would-be classical topic of algebraic 
geometry, which, however, did not attract much attention so far. We believe that time has come for its investigation 
and this paper is just a modest example of how one can approach the fundamental problems of this kind: very 
simple methods are quite effective and produce answers, which are not easy to foresee, and numbers, which are not 
easy to guess. This looks like a real and wonderful science to do. 



12 Appendix. Some elementary MAPLE programs for UMS studies 

We did our best to illustrate quantitative considerations of Universal Mandelbrot Set and its particular sections 
with modest illustrations. However, the number of illustrations in a printed text is necessarily restricted and can 
be non-sufficient for full visualization of the object. In order to cure this problem we collect in this appendix a set 
of sample MAPLE programs, which were used to generate some illustrations in the text. One can easily play with 
these simple programs, change parameters, accuracy of calculation and output formats in order to extract more 
information, numerical and visual. Programs are super-primitive, transparent and easy to modify, they work fast 
and smoothly on ordinary PC's. One can straightforwardly copy them into MAPLE file (with .raws extension) 
and use. When substituting desired parameters instead of the question marks, one should better do it in rational 
rather than decimal form, say a — 1/10 rather than a = 0.1. 

12.1 Cardioids 

MAPLE program for cardioid studies consists of just four lines: 
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> r:=l: 

> a:=?: b:=-(l+2*a)/3; 

> f :=r*(exp(I*t)+a*exp(I*2*t)+b*exp(I*3*t)) ; 

> plot( [Re(f ) ,Im(f ) ,t=-Pi. .Pi] , scaling=CONSTRAINED) ; 

(cubic case is presented, generalization is obvious). It remains to substitute various a and b instead of the question 
marks (say, a — 0.4 + 0.2 * /) and enjoy the pictures. For looking at more details, especially at the critical values 
t = and t = 7r, where cusps can occur (or t — J^y in general case) one can enhance resolution: 

> plot ( [Re (c) ,Im(c) ,t=Pi-0.01 . .Pi+0.01] , scaling=constrained) ; 

> plot( [t ,Re(c)/Im(c) ,t=Pi-0.01 . .Pi+0.01] ) ; 

> plot( [t,Im(c)/Re(c) ,t=-0.01. .+0.01]) ; 

Examples of output of this program are shown in Figs [5] and [6l 

12.2 UMS through discriminants and resultants 

> Fl:=f(x)-x: 

> F2:=f (f (x))-x: 

> F3:=f (f (f (x)))-x: 

> F4:=f (f (f (f (x))))-x: 

> F5:=f (f (f (f (f (x)))))-x: 

> F6:=f (f (f (f (f (f (x))))))-x: 

> G1:=F1: 

> G2:=simplify(F2/Gl) : 

> G3:=simplify(F3/Gl) : 

> G4:=simplify(F4/(G2*Gl))) ; 

> G6:=simplify(F6/(Gl*G2*G3)) ; 

> . . . 

> D2:=discrim(G2,x) ; 

> D3:=discrim(G3,x) ; 

> . . . 

> R24 : =resultant (G2 , G4 , x) ; 

> R36:=resultant(G3,G6,x) ; 

> ... 

12.3 Domains (1), (2) and (2,1) of M ax t +{1 ^ a)x 2 +c 

Parameter M in the program defines the number T of points in the picture. The bigger M the more detailed will 
be the plot, but computer time will also increase. To make sure that the program is working we added the line 
print( v k = ", k), one can safely omit it. 

> with(plots) : 
> 

> unassignCa' , 'b' ,'z' ,'t' ,'c') : 
> 

> a:=?: 

> b:=l-a: 
> 

> Dl : =f actor (discrim(a*x~3+b*x~2+c-x,x) ) ; 

> R21 :=f actor (resultant (a~3*x~6+2*a~2*x~5*b+a~2*x~4+2*a~2*x"3*c+a*b~2*x~4+2*a*x~3*b+a*x~2+2*x~2*c*a*b+x*c*a+ 

> D2:=f actor (simplify (discrim(a"3*x~6+2*a"2*x~5*b+a"2*x~4+2*a~2*x~3*c+a*b~2*x~4+2*a*x~3*b+a*x~2+2*x"2*c*a*b+ 
> 

> ## various choices of s and MID 

> #s : =evalf (solve (Dl , c) ) : 

> s : =evalf (solve (D2 , c) ) ; 

> #s:=evalf (solve (R21,c)) ; 

> MID:=s[l] ; 

> #MID:=(s[l]+s[2])/2. ; 
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> 

> P:= x -> a*x~3+b*x~2: 

> u:=exp(I*t) : 
> 

> zp:=(x,t)->(-b+root [2] (b~2+3*a*exp(I*t)/(3*a*x~2+2*b*x)))/(3*a) : 

> zm: = (x,t)->(-b-root [2] (b~2+3*a*exp(I*t)/(3*a*x~2+2*b*x)))/(3*a) : 

> sp:=solve(P(zp(x,t))+zp(x,t)-P(x)-x,x) : 

> sm:=solve(P(zm(x,t))+zm(x,t)-P(x)-x,x) : 
> 

> Tp:=0: Tm:=0: 

> M:=200: 

> for k to M do 
> 

> t:=evalf (2*Pi*k/M) : 

> N : =ArrayNumElems (Array ( [sp] ) ) ; 

> for i to N do 
> 

> wp:=allvalues(sp[i] ) : wm:=allvalues(sm[i] ) : 

> n:=ArrayNumElems (Array ( [wp] )) : # nm:=ArrayNumElems (Array ( [wm] )) : print (n,nm) ; 

> for j to n do 
> 

> Tp:=Tp+l: Tm:=Tm+l: 

> if n >1 then 

> Xp:=evalf (wp [j] ) : Xm: =evalf (wm [j] ) : 

> else 

> Xp : =evalf (wp) : Xm:=evalf (wm) : 

> end if : 

> Pp[Tp] :=evalf (zp(Xp,t)-(a*Xp~3+b*Xp~2)) : 

> Pm[Tm] :=evalf (zm(Xm,t)-(a*Xm"3+b*Xm~2)) : 
> 

> xpl:=Xp: xml:=Xm: 

> zpl:=a*xpl~3+b*xpl~2+Pp[Tp] : zml :=a*xml~3+b*xml~2+Pm[Tm] : 

> chp:=a*zpl~3+b*zpl~2+Pp[Tp] -xpl : 

> dm:=a*zml~3+b*zml~2+Pm[Tm] -xml : 

> ap:= evalf (Re (chp) ~2+Im(chp) ~2) : am: =evalf (Re(chm) ~2+Im(chm) ~2) : 
> 

> # MAGNIFY (Enhanced resolution for vicinity of a chosen value of 'c') 

> ## CENTER POSITION 

> zz:=s [1] ; 

> ### version of defining zz 

> #zz:=MID+I*0. : 

> ## RADIUS 

> rr:=0.3; 

> if rr>0 then 

> if (ap>10~(-5)) or abs (Pp [Tp] -zz) >rr then 

> Tp:=Tp-l: 

> else 

> fi: 

> if (am>10"(-5)) or abs(Pm[Tm] -zz)>rr then 

> Tm:=Tm-l: 

> else 

> fi: 

> else 

> if (ap>10~(-5)) then Tp:=Tp-l: fi: 

> if (am>10"(-5)) then Tm:=Tm-l: fi: 

> fi: 
> 

> od: 
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> od: 

> od: 
> 

> pp:=pointplot({seq([Re(Pp[n]) , Im(Pp [n] )] ,n=l. .Tp)}, 

scaling=CONSTRAINED , color=red, symbol=circle , symbols ize=5) : 

> pm:=pointplot({seq([Re(Pm[n]) , Im(Pm[n] )] ,n=l. .Tm)}, 

scaling=CONSTRAINED, color=red, symbol=circle, symbolsize=5) : 

> 

> display ({pp}, {pm}) ; a; 

12.4 3D tubes 

> unassignCa' , 'b' , 'u' ,'z' ,'t' , 'c') : 
> 

> a:=b->b~3: 

> c:=b->l. : 

> # I f ' I VALUE 

> MD:=1: 

> sp:=solve(4*a(b)*x~3+3*b*x~2+2*c(b)*x-MD*exp(I*t) ,x) : 
> 

> Tp:=0: Tm:=0: 

> M:=00: 

> M1:=15:M2:=60: 

> zmi := .2: zma : = . 8 : 

> for kl to Ml do 

> print ("k=" ,kl) ; 

> for k2 to M2 do 
> 

> t:=evalf (2*Pi*kl/Ml) : 

> b:=zmi+(zma-zmi)*k2/M2: 
> 

> N : =ArrayNumElems (Array ( [sp] ) ) ; 

> for i to N do 
> 

> wp : =allvalues (sp [i] ) : 

> n : =ArrayNumElems (Array ( [wp] ) ) : 

> for j to n do 

> Tp:=Tp+l: 

> if n >1 then 

> Xp:=evalf (wp[j] ) : 

> else 

> Xp:=evalf (wp) : 

> end if : 

> u:=evalf (Xp-(a(b)*Xp-4+b*Xp~3+c(b)*Xp~2)) : 

> Pp [Tp] : =arr ay ( [Re (u) , Im (u) , b] ) : 

> xpl:=Xp: 

> zpl:=a(b)*xpl~4+b*xpl~3+c(b)*xpl-2+Pp[Tp] [l]+I*Pp[Tp] [2] : 

> chp:=a(b)*xpl~4+b*xpl~3+c(b)*xpl~2+Pp[Tp] [l]+I*Pp[Tp] [2]-xpl: 

> ap:= evalf (Re (chp) ~2+Im(chp) ~2) : 

> if (ap>10~(-5)) then 

> Tp:=Tp-l: 

> fi: 

> od: 

> od: 

> od: 

> od: 
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> # PREPARE ARRAY FDR 3D PLOT 

> L:=l: 

> N:=Tp+Tm; 

> B:=array(l . . M) : 

> k:=0:j:=0: 
> 

> for i to Tp do 

> k:=k+l: 

> B[k] :=Pp[i] ; 

> od: 
> 

> for i to Tm do 

> k:=k+l: 

> B [k] : =Pm [i] ; 

> od: 
> 

> j:=j+l: 

> print(k,j,B[k]); 



> # PLOT 

> with(linalg) : 

> with(plots) : 

> with(plottools) : 

> setoptions3d(color=BLUE,symbol=CR0SS,symbolsize=3) ; 

> p:=pointplot3d(B,axes=B0XED) : 

> display (p) ; 



12.5 Fragments of Julia sheaf J a x' i +(i-a)x 2 +c- orbits of orders 1 and 2 vs c and a 

> unassign('x' , 'a' , 'b' , 'c') : 

> f :=x->a*x~3+b*x~2+c; 

> simplify (diff (f(f(x)) ,x)) ; 

> fpl:=x->diff (f (x) ,x) : 

> fp:=x->diff (f (f (x)) ,x) : 

> Gl:=f(x)-x; 

> F2:=f (f (x))-x: 

> G2:=simplify(F2/Gl) ; 



> ################################################ 

> a:=l/10; 

> b:=l.-a; 

> 

> Dl : =f actor (discrim(a*x~3+b*x~2+c-x,x) ) : 

> R21 :=f actor (resultant (a"3*x~6+2*a"2*x~5*b+a"2*x"4+2*a"2*x"3*c+ 
a*b~2*x~4+2*a*x~3*b+a*x~2+2*x~2*c*a*b+x*c*a+a*c~2+b~2*x~2+x*b+c*b+l,c-x+a*x~3+b*x~2,x)) : 

> D2:=f actor (discrim(a~3*x~6+2*a~2*x~5*b+a~2*x~4+2*a~2*x~3*c+ 
a*b~2*x~4+2*a*x"3*b+a*x"2+2*x"2*c*a*b+x*c*a+a*c"2+b~2*x~2+x*b+c*b+l,x))/R21 : 

> 

> # GET MIDDLE POINT 

> s : =evalf (solve (Dl , c) ) : 

> ## versions of defining 's' 

> #s:=evalf (solve(D2,c)) ; 
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> #s:=evalf (solve (R21,c)) ; 

> MID: = (s[l]+s[2])/2. ; 
> 

> # CHOOSE C VALUE 

> c: =-6. 24+1*0. : 

> ## version of defining 'c' 

> #c:=MID+I*0. ; 
> 

> sl:=solve(Gl,x) ; 

> s2:=solve(G2,x) ; 
> 

> Nl : =ArrayNumElems (Array ( [si] ) ) ; 

> N:=ArrayNumElems (Array ( [s2] )) ; 
> 

> # GET PAIRS 

> k:=0: 

> for i to N do 

> for j from i+1 to N do 

> if abs(f (s2[i])-s2[j])<0.0001 then 

> k:=k+l: 

> P[k] [1] :=i: 

> P[k] [2] :=j: 

> fi: 

> od: 

> od: 
> 

> print (P) ; 



> # SHOW ROOT POSITION 

> with(plots) : 

> pO:=pointplot({[Re(sl[l]) ,Im(sl[l])] , [Re (si [2] ) , Im(sl [2] )] , [Re (si [3] ) , Im(sl [3] )] }, 
color=BLACK, symbol=CR0SS , symbols ize= 15) : 

> pll:=pointplot({[Re(s2[P[l] [1] ] ) , Im(s2 [P [1] [1] ] )] } , color=red) : 

> pl2:=pointplot({[Re(s2[P[l] [2] ] ) , Im(s2 [P [1] [2] ] )] } , color=red) : 

> p21 : =pointplot ({ [Re (s2 [P [2] [1] ] ) , Im(s2 [P [2] [1] ] ) ] } , color=green) : 

> p22:=pointplot({[Re(s2[P[2] [2] ] ) , Im(s2 [P [2] [2] ] )] } , color=green) : 

> p31:=pointplot({[Re(s2[P[3] [1] ] ) , Im(s2 [P [3] [1] ] )] } , color=blue) : 

> P 32:=pointplot({[Re(s2[P[3] [2] ] ) , Im(s2 [P [3] [2] ] )] } , color=blue) : 

> display({p0,pll,pl2,p21,p22,p31,p32}) ; 

12.5.1 Stability of orbits 

> # GET STABILITY INFO 

> print ("ORDER 1") ; 

> Fl:=fpl(x): 

> for i to N/2 do 

> x:=sl [i] ; 

> print(abs(Fl) ,x) ; 

> od: 

> unassign( 'x' ) ; 
> 

> print ("ORDER 2") ; 

> F:=fp(x): 

> ## versions of defining 'F' 

> #F:=2*b*x; 

> #F:=3*a*x~2+2*b*x; 



62 



> 

> for i to N/2 do 

> x:=s2[P[i] [1]] ; 

> print (abs (F) ,x,f(x)) ; 

> od: 

> imassign( 'x' ) ; 



12.5.2 Attraction pattern 

> unassign('a' , 'b' , 'c') : 

> f :=x->a*x~3+b*x~2+c; 

> F2:=factor(f (f (x))-x) ; 
> 

> a:=l/3: b:=l-a: 

> R21 :=f actor (resultant (a~3*x~6+2*a~2*x~5*b+a~2*x~4+2*a~2*x~3*c+ 
a*b"2*x"4+2*a*x"3*b+a*x~2+2*x~2*c*a*b+x*c*a+a*c~2+b~2*x~2+x*b+c*b+l,c-x+a*x~3+b*x~2,x)) ; 

> D2:=f actor (discrim(a~3*x~6+2*a~2*x~5*b+a~2*x~4+2*a~2*x~3*c+ 
a*b"2*x"4+2*a*x~3*b+a*x~2+2*x~2*c*a*b+x*c*a+a*c~2+b~2*x~2+x*b+c*b+l,x)) ; 

> 

> rl2 : =evalf (solve (R21 , c) ) ; 

> d2:=evalf (solve (D2,c)) ; 



> ND:=4: 

> k:=0: 

> for t to ND+1 do 

> c:=d2[3]+0.01*exp(I*2*Pi/ND*(t-l)+I*Pi/2) ; 

> C[t] :=c; 

> sx : =evalf (solve (a~3*x~6+2*b*a~2*x~5+a~2*x~4+2*a~2*x~3*c+a*b~2*x~4+ 
2*b*a*x"3+a*x"2+2*b*a*x~2*c+a*x*c+a*c~2+b~2*x~2+b*x+c*b+l,x)) ; 

> for j to 6 do 

> R[(t-l)*6+j] :=sx[j] ; 

> k:=k+l: 

> od: 

> print ("t=" ,t) ; 

> od: 

> print (k) ; 

> # COMPUTE PATHS WITH RANDOM START 

> rf :=rand(-100. .100) : 

> L:=l: 

> c:=evalf (C[L+1]) ; 

> k:=0: 

> BL:=1:EL:=20: 

> LM:=2.8: 

> for i from BL to EL do 
> 

> cx:=rf ()/30+I*rf O/30: 

> CX:=cx: 

> N:=l: 

> jl:=0: 

> ep: =0.0003: 

> for j to N do 

> cx:=evalf (CX+0*exp(I*2*Pi*j/N)*ep) : 

> abs(cx-CX) ; 

> for il to 50 do 

> cx : =f (cx) ; 
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> if abs(Re(cx))<LM and abs(Im(cx))<LM then 

> k:=k+l; 

> S [k] :=cx; 

> fi: 

> od: 

> od: 
> 

> od: 

> print (k) ; 



> # PLOT DATA 

> with(plots) : 

> pr :=pointplot ({seq( [Re(R[n] ) ,Im(R[n] )] ,n=L*6+l . .L*6+6)>, 
scaling=CONSTRAINED , color=blue , symbol=cross , symbolsize=25) : 

> psl:=pointplot({seq([Re(S[n]) ,Im(S[n])] ,n=l. .k)}, 
scaling=CONSTRAINED , color=red , symbol=circle , symbols ize=5) : 

> display ({pr,psl>) ; 
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